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Abstract 


This  technical  report  presents  theoretical  and  experimental  results  of  phase  I  of  the 
Technology  Investment  Fund  (TIF)  project  on  low  probability  of  intercept,  that  is 
bistatic  SAR.  In  bistatic  SAR,  the  transmitter  and  receiver  are  spatially  separated 
and  hence,  the  risks  of  detection  and  localization  are  significantly  reduced,  i.e.,  its 
vulnerability  to  jamming  is  reduced  and  its  survivability  significantly  increased. 
Bistatic  SAR  images  include  information  that  is  complementary  to  monostatic  im¬ 
ages  due  to  the  different  scattering  mechanisms  involved.  Such  information  could 
lead  to  the  development  of  new  techniques  for  automatic  target  recognition  and 
classification. 

The  project  investigates  the  feasibility  of  bistatic  SAR  and  identifies  performance 
limits  through  a  trade-off  analysis  between  radar  parameters/geometry  and  achiev¬ 
able  resolution.  As  the  main  result,  conclusions  are  drawn  regarding  the  bistatic 
observation  time  and  the  subsequent  imaging  performance  under  different  bistatic 
configurations,  as  well  as  the  performance  degradation  due  to  severe  oscillator 
phase  noise  or  jitter.  The  report  also  includes  a  description  of  a  bistatic  SAR  sim¬ 
ulator,  an  implemented  version  of  a  time  domain  bistatic  SAR  processor,  and  the 
analysis  of  existing  experimental  bistatic  clutter  data  sets. 


Resume 


Le  present  rapport  technique  donne  les  resultats  theoriques  et  experimentaux  de  la 
phase  I  d’un  projet  mene  dans  le  cadre  du  Fonds  d’investissement  technologique 
(FIT)  sur  la  faible  probability  d’ interception  (FPI),  c’est-a-dire  le  radar  a  synthese 
d’ouverture  (RSO)  bistatique.  Dans  le  RSO  bistatique,  l’emetteur  et  le  recepteur 
sont  separes  spatialement,  ce  qui  reduit  significativement  les  risques  de  detection 
et  de  localisation,  c’est-a-dire  que  la  vulnerability  du  radar  au  brouillage  est  reduite 
et  sa  surviabilite,  fortement  accrue.  Les  images  de  RSO  bistatique  component  de 
F information  qui  vient  completer  les  images  de  radar  monostatique  en  raison  des 
differents  mecanismes  de  diffusion  en  jeu.  Cette  information  pourrait  mener  a  la 
mise  au  point  de  nouvelles  techniques  de  reconnaissance  et  de  classification  auto- 
matiques  des  cibles. 

Le  projet  comporte  une  etude  de  faisabilite  du  RSO  bistatique  et  en  definit  les  li- 
mites  de  rendement  par  une  analyse  comparative  entre  les  parametres  et  la  geometric 
du  radar,  d’une  part,  et  la  resolution  realisable.  Comme  grand  resultat,  des  conclu¬ 
sions  sont  tirees  au  sujet  du  temps  d’observation  du  radar  bistatique  et  l’obtention 
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subsequente  d’ images  en  fonction  de  differentes  configurations  de  radar  bistatique. 
ainsi  qu’au  sujet  de  la  degradation  du  rendement  a  cause  de  grave  bruit  de  phase 
d’oscillateur  ou  gigue.  Le  rapport  comprend  aussi  une  description  d’un  simulateur 
de  RSO  bistatique,  une  version  mise  en  aeuvre  d'un  processeur  de  RSO  bistatique 
temporel  et  Panalyse  des  series  de  donnees  disponibles  de  clutter  de  radar  bistatique 
experimental. 


DRDC  Ottawa  TR  2004-190 


Executive  summary 


In  military  environments,  imaging  based  on  monostatic  synthetic  aperture  radar 
(SAR)  is,  due  to  its  active  illumination  of  the  scene,  easily  detected  and  thus  highly 
vulnerable  to  electronic  counter  measures  or  even  destruction.  Use  of  a  bistatic  (or 
multi-static)  system,  with  its  transmitter  spatially  separated  from  the  ‘silent’  receiv¬ 
ing  systems  significantly  reduces  the  risk  of  receiver  detection  and  increases  surviv¬ 
ability.  In  addition  to  low  probability  of  intercept  (LP1),  such  a  bistatic  SAR  config¬ 
uration  offers  additional  surveillance  capabilities,  such  as  complementary  imagery 
due  to  forward-scattering  measurements.  The  proposed  concept  opens  a  wide  field 
of  future,  related  research  and  will  even  lead  to  different  applications,  like  bistatic 
SAR  in  combination  with  ground  moving  target  indication  or  topographic  mapping. 
Even  forward  looking  SAR  becomes  possible  with  bistatic  configurations  (impos¬ 
sible  with  monostatic  systems)  which  might  be  exploited  in  landing  aids  for  aircraft 
or  obstruction  warning  devices  in  helicopters.  This  work  could  also  serve  as  foun¬ 
dation  of  R&D  in  the  field  of  covert  radar  detection  and  passive  imaging,  which  also 
exploits  bistatic  geometry  but  uses  a  non-cooperative  illuminator  of  opportunity. 

The  purpose  of  this  Technology  Investment  Fund  (TIF)  project  is  to  investigate  the 
feasibility  and  the  limitations  of  a  bistatic  SAR  system.  This  includes  theoretical 
study  and  simulation  of  optimum  bistatic  SAR  configurations,  the  development  of 
robust  processing  algorithms  for  imaging,  and  an  experimental  demonstration  with 
the  upgraded  airborne  SpotSAR  system  and  a  specially  constructed  stationary  re¬ 
ceiver. 

This  report  presents  a  theoretical  framework  for  performance  analysis  prediction 
for  any  bistatic  geometry  with  respect  to  the  geometric  and  radiometric  ground  res¬ 
olution.  It  introduces  analysis  of  the  most  promising  configurations,  including  their 
performance  limits  and  relevant  problems,  as  well  as  the  development  and  imple¬ 
mentation  of  a  bistatic  SAR  simulator,  which  will  be  used  to  confirm  the  analysis 
and  to  test  bistatic  SAR  processors.  Finally,  the  report  describes  the  successful 
coherent  processing  of  two  existing  bistatic  clutter  data  sets  (space-borne  and  air¬ 
borne  data)  measured  at  the  Air  Force  Research  Laboratories  (AFRL)  in  Rome,  NY 
in  2002. 

This  research  provides  a  deep  level  of  understanding  of  the  capabilities  and  difficul¬ 
ties  associated  with  bi-  and  multistatic  SAR  concepts.  This  knowledge  will,  for  in¬ 
stance,  be  essential  to  perform  system  design  studies  and  accurate  cost/performance 
trade-off  analyses  of  SBR  designs  in  combination  with  airborne  systems  for  future 
military  SAR  applications.  The  application  of  such  concepts  could  significantly 
enhance  Canada’s  surveillance  capabilities  in  times  of  conflict  by  reducing  vulner¬ 
ability  to  receiver  jamming  and  to  direct  attack  of  the  surveillance  asset. 


iii 
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Presently,  Canada  is  in  the  favourable  situation  of  being  able  to  collect  data  from  na¬ 
tional  assets  including  both  space-borne  sensors,  RADARSAT-1  and  soon  RADAR- 
SAT-2,  and  experimental  airborne  SAR  systems,  such  as  DRDC’s  SpotSAR  or  En¬ 
vironment  Canada’s  C-band  SAR.  This  rather  unique  position  offers  the  opportu¬ 
nity  to  test  sophisticated  new  radar  concepts  using  measured  data  gathered  under 
realistic  military  conditions. 


Christoph  H.  Gierull;  2004;  Bistatic  Synthetic  Aperture  Radar;  DRDC  Ottawa 
TR  2004-190;  Defence  R&D  Canada  -  Ottawa. 
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Sommaire 


Dans  les  milieux  militaires,  l’imagerie  fondee  sur  le  radar  a  synthese  d’ouver- 
ture  (RSO)  monostatique  est,  en  raison  de  l’illumination  active  de  la  scene,  fa- 
cilement  detectee  et,  par  consequent,  hautement  vulnerable  aux  contre-mesures 
electroniques  ou  meme  a  la  destruction.  L’utilisation  d’un  systeme  bistatique  (ou 
multistatique),  dont  l’emetteur  est  spatialement  separe  des  systemes  de  reception 
‘silencieux’,  reduit  significativement  le  risque  de  detection  du  recepteur  et  accroit 
la  surviabilite.  Outre  la  faible  probability  d’ interception  (FPI),  une  telle  configura¬ 
tion  de  RSO  bistatique  offre  des  capacites  additionnelles  de  surveillance,  corame 
1’imagerie  complementaire  grace  a  des  mesures  de  prodiffusion.  Le  concept  pro¬ 
pose  ouvre  un  vaste  domaine  de  recherches  connexes  et  menera  meme  a  differentes 
applications,  comme  l’utilisation  du  RSO  bistatique  en  conjonction  avec  la  car¬ 
tographic  topographique  ou  l’indication  de  cibles  mobiles  au  sol.  Meme  le  RSO 
a  visualisation  vers  l’avant  devient  possible  avec  les  configurations  de  radar  bi¬ 
statique  (impossible  avec  les  systemes  monostatiques),  et  il  pourrait  etre  exploite 
comme  aide  a  l’atterrissage  pour  les  dispositifs  de  detection  d’obstacle  ou  d’aeronef 
a  bord  d’helicopteres.  Ce  travail  pourrait  aussi  servir  de  fondement  a  la  recherche 
et  developpement  dans  le  domaine  de  la  detection  radar  discrete  et  de  l’imagerie 
passive,  qui  exploite  aussi  la  geometrie  bistatique,  mais  qui  utilise  un  illuminateur 
de  choix  non  cooperatif. 

Le  but  du  present  projet  mene  dans  le  cadre  du  Fonds  d’investissement  technolo- 
gique  (FIT)  est  l’etude  de  la  faisabilite  et  des  limites  d’un  systeme  RSO  bistatique, 
y  compris  la  simulation  et  l’etude  theorique  de  configurations  optimales  de  RSO 
bistatique,  la  mise  au  point  de  solides  algorithmes  de  traitement  d’imagerie  et  une 
demonstration  experimental  au  moyen  du  systeme  SpotSAR  aeroporte  ameliore  et 
d’un  recepteur  stationnaire  particular. 

Le  present  rapport  donne  un  cadre  theorique  de  predication  d’ analyse  du  rendement 
de  n’importe  quelle  geometrie  bistatique  en  ce  qui  conceme  la  resolution  au  sol 
geometrique  et  radiometrique.  II  introduit  l'analyse  des  configurations  les  plus  pro- 
metteuses,  y  compris  leurs  limites  de  rendement  et  les  problemes  pertinents,  ainsi 
que  le  developpement  et  la  mise  en  place  d’un  simulateur  de  RSO  bistatique,  qui 
servira  a  confirmer  l'analyse  et  a  mettre  a  l’essai  les  processeurs  de  RSO  bistatique. 
Enfin,  le  rapport  decrit  le  traitement  coherent  reussi  de  deux  series  disponibles  de 
donnees  de  clutter  bistatique  (donnees  spatiales  et  de  vol)  prises  a  l’Air  Force  Re¬ 
search  Laboratory  (AFRL)  a  Rome  (N.Y.)  en  2002. 

Cette  recherche  permet  d’obtenir  une  comprehension  en  profondeur  des  capacites 
et  des  difficultes  associees  aux  concepts  du  ROS  bistatique  et  multistatique.  Ces 
connaissances  sont,  par  exemple,  essentielles  a  la  conduite  d’etudes  de  conception 
de  systemes  et  d’analyses  comparatives  cots/rendement  precises  des  conceptions 
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du  radar  satellise  (RS)  en  conjonction  avec  les  systemes  aeroportes  en  vue  dupli¬ 
cations  militaires  futures  du  RSO.  L’ application  de  ces  concepts  pourrait  ameliorer 
significativement  les  capacites  de  surveillance  du  Canada  en  periodes  de  conflit.  en 
reduisant  la  vulnerability  au  brouillage  des  recepteurs  et  aux  attaques  directes  des 
dispositifs  de  surveillance. 

Le  Canada  se  trouve  actuellement  dans  une  situation  favorable,  ctant  en  mesure 
de  recueillir  des  donnces  a  partir  de  dispositifs  nationaux,  y  compris  les  deux  cap- 
teurs  spatiaux  (RADARSAT-1  et  bientt  RADARSAT-2)  et  des  systemes  spatiaux 
experimentaux  (comme  SporSAR  dc  RDDC  ou  le  RSO  en  bande  C  d'Environnc- 
ment  Canada).  Cette  position  plutt  unique  offre  la  possibility  de  mettre  a  l’essai 
de  nouveaux  concepts  perfectionnes  dans  le  domaine  du  radar  a  l’aide  de  donnees 
mesurees  recueillies  dans  des  conditions  militaires  rcalistes. 

Christoph  H.  Gierull;  2004;  Radar  a  synthese  d’ouverture  bistatique;  DRDC  Ottawa 
TR  2004-190;  R  &  D  pour  la  defense  Canada  -  Ottawa. 
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1  Introduction 


In  military  environments,  monostatic  synthetic  aperture  radar  (SAR)  is,  due  to  its 
active  illumination  of  the  scene,  easily  detected  and  thus  highly  vulnerable  to  elec¬ 
tronic  counter  measures  or  even  destruction.  Use  of  a  bistatic/multi-static  system, 
with  its  transmitter  spatially  separated  from  the  ‘silent’  receiving  systems  signifi¬ 
cantly  reduces  the  risk  of  receiver  detection  and  increases  survivability.  Potentially 
qualified  candidates  for  the  transmitter  are  space  based  radars  (SBR)  or  aircrafts 
flying  at  high  altitudes.  The  receivers  can  be  mounted  on  highly  maneuverable  air¬ 
craft  or  UAVs  that  may  image  the  scene  at  different  times  and  from  different  angles. 
In  addition  to  low  probability  of  detection,  such  a  bistatic/multistatic  SAR  config¬ 
uration  offers  additional  surveillance  capabilities,  such  as  complementary  imagery 
due  to  forward-scattering  measurements.  This  may  lead  to  the  development  of  new 
techniques  for  automatic  target  recognition  and  classification  and  could  enhance 
the  visibility  of  stealthy  targets.  Furthermore,  improved  ground  moving  target  in¬ 
dication  (GMTI)  and  topographical  surveying  and  mapping  of  the  scene  become 
possible. 

Recently,  there  has  been  a  growing  interest  in  the  use  of  space-based  radars  (SBR) 
to  support  the  future  military  requirements  of  Canada  and  its  allies.  These  require¬ 
ments  include  24/7  wide-area  surveillance  of  land  targets  under  all  weather  condi¬ 
tions.  Consequently,  future  military  SBR  systems  will  be  equipped  with  synthetic 
aperture  radar.  Wide  area  coverage  from  space  by  SAR  requires  resolution  to  be 
traded  off  with  swath  coverage.  It  is  further  limited  by  the  antenna  power  handling 
capacity,  limited  down-link  capacity,  lower  SNR  (due  to  larger  distances  between 
sensor  and  targets)  and  a  less  flexible  imaging  geometry  compared  to  airborne  SAR 
systems.  On  the  other  hand,  airborne  SAR  is  highly  vulnerable  because,  as  an  active 
system,  its  transmitted  energy  can  be  detected  and  hence  intercepted,  particularly 
when  operating  in  a  hostile  environment.  Interception  includes  active  receiver  jam¬ 
ming  as  well  as  the  actual  destruction  of  the  radar  platform. 

This  vulnerability  may  be  overcome  during  a  conflict  by  using  a  bistatic  ‘low  prob¬ 
ability  of  intercept’  SAR  concept,  where  the  radar  transmitter  may  be  located  on 
a  satellite,  while  the  receiver  is  located  on  an  airborne  platform,  such  as  a  UAV. 
Since  the  receiving  radar  system  operates  in  a  passive  (‘silent’)  mode,  the  risk  of 
receiver  detection  and  localization  is  significantly  reduced,  i.e.,  its  vulnerability 
to  jamming  is  reduced  and  its  survivability  significantly  increased.  Compared  to 
monostatic  space-based  SAR,  a  higher  SNR  of  the  scattered  echoes  can  also  be 
achieved  because  of  the  smaller  distances  between  the  targets  and  the  receiver.  Fur¬ 
thermore,  complementary  information  in  the  SAR  images  can  be  expected  due  to 
the  forward-scattering  measurements  of  bistatic  systems.  For  instance,  the  roofs 
of  buildings  might  be  the  major  source  of  scattering  in  bistatic  SAR  imagery  com- 
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pared  to  the  walls  in  conventional,  monostatic  SAR  images.  Combined  with  the 
advantages  of  a  more  flexible  imaging  geometry,  it  can  be  seen  how  such  a  bistatic 
SAR  concept  can  enhance  surveillance  capabilities  such  as  ground  moving  target 
indication  (GMTI)  and  topographical  surveying  and  mapping. 

In  contrast  to  conventional,  monostatic  SAR.  the  required  bistatic  configuration  has 
a  far  more  complex  geometry  that  leads  to  certain  unique  problems  that  must  be 
dealt  with  during  image  formation.  Examples  of  such  problems  include  greater 
variation  of  Doppler  frequency  with  range,  non-linear  Doppler  frequency  modula¬ 
tion  in  the  along-track  direction,  reception  of  transmitter  antenna  mainlobe  energy 
through  receiver  antenna  sidelobes,  and  fold-over  effects  resulting  from  range  am¬ 
biguities.  In  order  to  generate  high-resolution  SAR  images,  the  bistatic  imaging 
geometry  must  be  known  very  accurately.  This  makes  great  demands  on  both  the 
motion  measurement  system/equipment,  such  as  INS  and  GPS,  and  the  SAR  signal 
processing/autofocus  algorithms.  Furthermore,  the  echo  pulses  must  be  coherently 
measured,  i.e.,  the  phase  information  of  the  transmitted  pulse  has  to  be  preserved. 
Hence,  the  problem  of  time  synchronization  between  the  two  systems  must  be  re¬ 
solved  because  the  transmitter  and  receiver  are  spatially  separated. 

Although  the  investigation  and  application  of  bistatic  radar  is  by  no  means  new 
[1,  2,  3,  4],  there  appears  to  be  a  resurgence  caused  by  technological  progress 
and  military  exigency,  particularly  for  imaging  radars.  Based  on  a  growing  con¬ 
fidence  of  technological  feasibility,  economical  affordability,  and  due  to  palpable 
advantages  of  bistatic  SAR  for  military  applications  worldwide,  scientific  activities 
have  intensified  recently.  Until  two  years  ago  only  a  handful  of  papers  were  pub¬ 
lished  in  the  open  literature.  At  the  international  conference  on  Geoscience  and 
Remote  SenSing  (IGARSS’03),  several  papers  about  bistatic  SAR  appeared.  One 
year  later,  at  the  European  SAR  Conference  EUSAR’04.  there  were  two  full  ses¬ 
sions  dedicated  to  phenomenology  and  processing  of  bistatic  and  multistatic  SAR 
data.  While  the  interest  in  bistatic  SAR  was  in  the  past  mostly  stamped  by  military 
advantages,  even  civil  applications  are  starting  to  be  discussed  [5,  6].  Furthermore, 
highly  sophisticated  bi-  and  multistatic  space-borne  system  concepts  are  being  in¬ 
vestigated  for  enhanced  applications  such  as  interferometry  [7,  8,  9]  and  SAR/MTI 
[10,  11,  12].  A  step  even  further  ahead  is  the  possibility  of  passive  (or  parasitic) 
SAR  imaging,  which  also  exploits  a  bistatic  geometry  but  using  a  non-cooperative 
illuminator  of  opportunity  [13,  14,  15]. 

However,  a  thorough  and  complete  theoretical  performance  analysis  with  regard 
to  the  maximal  geometric  resolution  determined  by  the  three-dimensional  flight 
geometry  and  the  maximum  observation  time  is  noticeably  absent  from  the  bistatic 
SAR  literature.  Further,  due  to  lack  of  necessity  in  monostatic  SAR  (with  only  one 
local  oscillator  (LO)  for  up-  and  down  conversion  of  the  signals)  the  radiometric 
deterioration  caused  by  the  severe  phase  noise  of  two  separate  LOs  has  not  been 
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extensively  investigated  in  the  past. 

This  TIF-project  is  divided  into  two  phases.  In  Phase  I,  the  theoretical  background 
of  bistatic  SAR  has  been  derived  via  stepwise  mathematical  modeling  of  the  re¬ 
ceived  echo  signals.  Simultaneously,  a  simulation  tool  and  bistatic  SAR  processing 
algorithms  are  developed  to  evaluate  and  demonstrate  the  performance  of  the  new 
SAR  processing  concepts.  A  detailed  investigation  of  bistatic  SAR  performance, 
like  achievable  geometric  and  radiometric  resolution,  has  been  performed  over  a 
wide  range  of  bistatic  geometries  and  radar  operating  parameters.  In  Phase  II  of  the 
project,  deeper  insight  into  the  functionality  of  a  bistatic  SAR  system  and  the  phys¬ 
ical  properties  of  its  SAR  imagery  will  be  gained  with  experimental  measurements 
conducted  by  DRDC  Ottawa’s  airborne  SpotSAR  radar  system. 
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2  Bistatic  SAR  signal  model 


In  this  Section  a  bistatic  SAR  signal  model  is  derived  as  well  as  the  predicted  SAR 
processing  degradation  due  to  oscillator  frequency  mismatch  between  the  two  sep¬ 
arated  transmit  and  receive  local  oscillators.  The  output  signal  from  two  ideally 
non-drifting  oscillators  would  be  a  pure  sine  wave  of  exactly  identical  frequency, 
but  any  real  device,  even  the  most  stable,  is  disturbed  by  unavoidable  processes 
such  as  random  noise  and  drifts  due  to  aging  and/or  environmental  effects. 

The  starting  point  for  the  analysis  of  bistatic  SAR  is  the  three-dimensional  geometry 
depicted  in  Fig.  1.  While  it  was  sufficient  for  investigation  of  bistatic  radar  in  the 
past  to  only  consider  two-dimensional  geometries,  analysis  of  bistatic  SAR  requires 
a  three-dimensional  geometry,  e.g.,  [16,  17].  As  a  convention  in  the  remainder  of 
this  document,  vector  quantities  are  underlined  and  matrices  are  written  in  bold. 
The  transpose  operator  is  denoted  as  superscript  complex  conjugate  as  *  and 
complex  conjugate  transpose  as  the  Hermitian  operator  f. 


Figure  1 :  Bistatic  SAR  geometry  with  constant  velocity  of  each  platform. 
Consider  two  independent,  spatially  separated,  radar  platforms  moving  along  straight 
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paths  in  the  three-dimensional  Cartesian  scene  co-ordinate  system,  which  is  de¬ 
noted  by  subscript  s,  appendix  D.l.  The  two  platforms  are  not  necessarily  required 
to  be  aircrafts.  For  example,  using  a  satellite  as  illuminator  is  a  particularly  inter¬ 
esting  configuration,  as  will  be  shown  later. 

A  detailed  definition  and  description  of  all  involved  co-ordinate  systems  is  given 
in  appendix  D.  The  center  (origin)  of  the  scene  is  defined  as  the  intersection  of 
the  transmitter  and  receiver  pointing  vectors  with  the  ground  plane  at  time  t  =  0. 
We  are  assuming  that  both  platforms  are  flying  in  a  plane  parallel  to  the  ground 
plane,  i.e.,  with  constant  altitudes  if,  and  with  constant  velocities  for  i  =  1,2. 
These  assumptions  are  not  restrictive  and  are  only  made  to  clarify  the  following 
derivations. 

The  range  vector  from  the  transmitter  (henceforth  denotes  by  subscript  1)  to  a  point 
r  on  the  ground  at  time  t  is  given  as 

Ei(t,r)  =  (^oi -r)  +VJ,  (1) 

where  Rq^  =  ^  (0, 0)  is  the  vector  from  the  scene  center  to  the  transmitter  at  time 
zero,  see  appendix  F.  This  vector  is  commonly  determined  by  the  depression  an¬ 
gle  in  elevation  and  the  squint  angle  in  azimuth  at  the  time  origin,  see  appendix 
A.  Analogously,  the  range  vector  for  the  receiver  is  denoted  by  subscript  2.  Note 
that  the  subscripts  1,2  are  used  rather  than,  e.g.,  ‘t,r’  to  indicate  the  complete  inter¬ 
changeability  of  receiver  and  transmitter,  from  a  physical  model  point-of-view. 

2.1  Constant  oscillator  frequency  offset 

2.1.1  Time  history 

At  subsequent  instances  in  time,  denoted  by  tn  =  nAT  for  n  =  0, 1, 2, . . .,  the 
transmitter  coherently  radiates  a  linear  frequency  modulated  (LFM)  pulse  with  car¬ 
rier  frequency  o>i  of  form  Re{strans(f)}  with 

Strans(f)  =  ^  (t_tn)V(u'lt+0l(o)+4l(t))  rect,  f-  l'n^\  +Tli(t),  (2) 

71=0  '  ' 

where  B  is  the  chirp  bandwidth,  T  is  the  pulse  duration,  and  AT  is  the  pulse  rep¬ 
etition  interval  (PRI).  The  term  ni(f)  represents  unavoidable  additive  noise,  e.g., 
the  thermal  sensor  noise  in  the  transmitter  hardware.  It  is  usually  modeled  as  a 
stationary  complex  normal  distributed  stochastic  process,  i.e.,  a  Gaussian  process 
with  expectation  zero  and  auto-covariance  function  Cnini  (t).  However,  since  we  are 
at  this  point  interested  in  phase  noise  effects,  the  additive  noise  term  is  omitted  in 
the  following  for  didactical  clarity.  The  stochastic  process  £i  (t)  denotes  the  ‘multi¬ 
plicative’  phase  noise,  i.e.,  describes  oscillator  phase  instabilities  such  as  jitter.  The 
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constant  initial  phase  ^i(O),  at  time  origin,  does  not  effect  the  SAR  performance 
and  can  be  set  to  zero. 

It  is  sufficient  for  the  range  processing  to  concentrate  on  one  particular  pulse  so  that 
the  summation  in  (2)  will  be  abandoned  henceforth  for  notational  simplicity.  The 
signal  is  reflected  from  a  scatterer  at  ground  location  r  and  impinges  on  the  receiver 
with  a  corresponding  time  delay  r: 

sref(f)  =  rect  ^  ~  ~  ,  (3) 

where  the  phase  noise  is  ignored  at  the  moment.  It  is  important  to  note  that,  al¬ 
though  the  time  delay  r  is  strictly  speaking  itself  time  dependent,  it  can  be  treated 
as  constant  over  the  duration  of  a  pulse  T  (i.e.,  r  =  r(£n))  but,  of  course,  will  vary 
slowly  from  pulse-to-pulse1. 

The  down-conversion  of  this  RF  signal  to  baseband  with  receiver  oscillator  fre¬ 
quency  U2  can  be  expressed  by  multiplication  with  exp(—  juj2t)  and  yields  the  re¬ 
ceived  signal 

srec(f)  =  e~juJ1T  (t-tn-rf  gjAu;t  rect  )  (4) 

where  A u  =  lji  —  u>2  denotes  the  constant  frequency  offset  between  the  two  inde¬ 
pendent  oscillators.  Note  that  a  constant  frequency  offset  means  a  linear  phase  drift. 
The  time  history  of  the  time-shifted  and  phase-rotated  received  pulses  is  depicted 
in  Fig.  2. 


PRI 


PRI 

- 


Figure  2 : 


k*i  +* 


Illustration  of  the  phase  and  time  shifts  of  the  received  pulses. 


t 


2.1.2  Fast-time  (range)  processing 

The  range-compression  is  performed  by  convolution  with  filter  function  h(t),  which, 
when  chosen  to  be  the  time-inverted  conjugate  of  the  demodulated  transmit  signal 

h(t)  =  exp  ^-j7r^f2^  rect  ,  (5) 

'This  may  be  exploited  to  achieve  high  resolution  in  the  along-track  or  Doppler  direction. 
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maximizes  the  signal-to-noise  ratio  (SNR).  The  matched  filter  output  or  impulse 
response  function  is  given  as 


OO  OC 

p(t)  =  (Srec  *  h)(t)  =  J  srec(t')h(t  -  t')dt'  =  J  srec(t  -  t')h{t')dt' 

—oo  —  OO 

OO 


B  ,12 


■e  j7TTl"  rect 


t  —  t'  —  tn  —  T 


rect 


dt'. 


(6) 


The  two  rectangular  functions  within  the  integrand  of  (6)  restrict  the  integration 
interval  to  the  time-dependent  values  a(t)  and  b(t)  (to  be  specified),  that  is: 


b(t) 

p[t)  =  e_?‘JlTejAa't  J  e J7rf((t-tn_'r)2_2(t_tn-r)t,)e-JAa,(' 


dt' 


a(t) 


—  g-JWlTgjAwtgj 


m 


,-j(2*§(t-tn-T)+Au)  t' 


o(t) 


b(Q  +  a(t) 


—  g-juiT^Auit^jnf  (t-tn-T)2e~j  (2ir§ (t-t„-r)+Au;) 

•  2(b(t)  —  a(f))  sine  ^(27t ^(f  -  t„  -  r)  +  Aw) 

The  definition  of  the  integration  interval  is  illustrated  in  Fig.  3. 


b(t)  —  a{t ) 


(7) 


Figure  3  :  Illustration  of  integration  interval  for  range-compression. 


The  integration  interval  b(t)  —  a(t)  is  non-zero  when  the  right  rectangle  starts  to 
touch  the  left  one  at  time  t  =  tn  +  r  +  T,  is  maximum  when  they  are  completely 
overlap  at  t  —  tn  +  r,  and  linearly  decreases  until  t  =  tn  +  t  —  T.  The  triangular 
time  history  can  be  expressed  analytically  as 

b(t)  -  a(t)  =  (T  -\t-tn-  t\)  rect  ^ ^  ^  .  (8) 
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The  sum  of  the  two  limits  can  be  expressed  as 


b(t )  +  a(t)  =  (t  —  tn  —  r)  rect 


t  tn  T 

2  T 


(9) 


From  (8)  it  is  clear  that  within  the  vicinity  of  the  reflected  pulses  t  =  tn  +  t  the  in¬ 
tegration  region  can  be  approximated  by  the  constant  pulse  length  T,  as  commonly 
done  in  the  literature: 


The  effect  of  an  oscillator  frequency  mismatch  becomes  apparent,  namely  the  posi¬ 
tion  of  the  peak  is  displaced  from  the  correct  location.  Setting  the  argument  of  the 
sinc-function  in  (10)  to  zero  yields 


to  =  tn  +  T - —  or  B(t0  -  tn  -  t)  =  TAf, 


(ID 


where  A /  =  Auj/2n  =  fi  —  fi-  As  the  range  resolution  is  determined  by  the 
inverse  of  the  pulse  bandwidth  (to  =  1  / B),  the  value  of  B(to  —  tn  —  r)  =  BAt 
represents  the  range  displacement  in  multiples  of  the  resolution  cell  size. 

2.1.3  Slow-time  (Doppler)  processing 

For  the  subsequent  Doppler  processing  of  a  target  at  ground  location  r,  the  pulse-to- 
pulse  time  history  is  analysed  along  the  corresponding  time  instances  tn  -I-  r(fn,  r). 
Selecting  the  values  of  the  impulse  response  function  (10)  at  these  times  yields  the 
phase  signal 


SD(tn,r)  =  exp  (  -  ju\T{tn,r))  exp  (jAu(tn  +  r(fn,r))) 

=  exp  (  -  j(u>i  -  Ao;)r(fn,r))  exp  (jAu>tn),  (12) 


where  the  magnitude  has  been  omitted  for  simplicity.  The  first  exponential  term 
shows  an  effective  carrier  frequency  shift  for  the  transmitter.  However,  this  shift 
is  negligible  because  reasonable  offsets  will  only  be  a  small  fraction  of  the  center 
frequency  |Au;/u;i|  <C  1.  In  contrast,  the  second  phase  term  cannot  be  neglected 
and  will  affect  the  Doppler  processing. 

By  omitting  the  offset,  i.e.  Aw  =  0,  the  first  exponential  term  in  (12)  represents 
the  common  phase  (Doppler)  history  caused  by  the  changing  distance  (time  delay) 
between  transmitter,  target  and  receiver  over  many  pulses 


exp  (  -  jvir(tn,r))  =  exp  (  -  jMIL8i(*n,r)||  +  lli^n, r)||)),  (13) 
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where  k\  =  2tt  f\/c  denotes  the  transmitted  wavenumber.  To  get  deeper  insight, 
it  is  sufficient  now  to  consider  a  quasi-monostatic  geometry  with  an  approximate 
quadratic  range-time  history.  It  will  be  shown  later  that  such  scenarios  comprise 
bistatic  flight  geometries  with  parallel  paths  of  both  platforms  as  well  as  either  a 
stationary  transmitter  or  receiver.  For  a  conventional  monostatic  SAR  with  identi¬ 
cal  distances,  equation  (13)  yields  the  classic  Doppler  phase  history  (second  order 
Taylor  series): 

Ml)  =  2&1||£(f,r)||  “  2h  ( \(f, )  +  -  *«)2)  .  (I4) 

where,  for  notational  simplicity,  the  slow-time  tn  has  been  replaced  by  t  and  k\  by 
k.  This  convention  will  be  used  henceforth  when  possible  without  confusion.  The 
time  where  the  platform  is  broadside  to  the  scatterer  located  at  r  is  denoted  by  4 
and  the  distance  between  them  (at  this  instant  in  time)  by  Rb.  For  non-squinted 
SAR  geometries  and  stationary  reflectors,  the  time  tb  coincides  with  the  point  of 
closest  approach  [18]. 


Inserting  (14)  into  (12)  yields  the  slow-time  or  Doppler  signal  that  can  be  written 
as 

/  Aw 

-j2k  lRb{r)  -  —tb 

exp  (• i2k  (t^  - h)  -  dfe(<  - 4)2))  ■  (15) 

where  constant  terms  will  be  omitted  from  here  onwards.  The  linear  coefficient 
is  tantamount  to  a  quasi  radial  velocity  vr  =  Au>/(2k).  Analogously  to  the  range- 
compression,  the  Doppler  signal  is  convolved  with  the  expected  signal  of  a  reflector 
at  location  4  =  0  (reference  function): 

OO 

Pd(4i)  =  (sD*sie{)(t)  =  J  sD(t  -  t')sre{(t')dt',  (16) 

—  OO 

where 

*„,((, r)  =  exp  ( 2)  rect  (£)  ■  07) 

where  Ts  denotes  the  synthetic  aperture  time  over  which  the  particular  reflector  is 
observed  by  the  radar.  The  maximum  duration  that  the  reflector  is  in  the  antenna 
beamwidth  is  usually  expressed  by  the  antenna  beamwidth  03<ib  and  the  slant  range 
distance  Ts  =  93dBRb(r)/va.  Using  this  value,  the  point  spread  function  reads 


so(t,r)  =  exp 


Po(t,r)  =  e 


=  ^  tb J  e  i2k ( ^ 


-T./2 


DRDC  Ottawa  TR  2004-190 


9 


Rb(r) 


(t  -  tb) 


(18) 


Again,  the  peak  is  displaced  by  A tb  =  (tb  —  tb),  where  t'b  is  computed  by  setting 
the  argument  of  the  sinc-function  to  zero 


tb  —  tb  + 


A  uRb(r) 
kv2a 


or  A tb  =  (tb  -  tb)  =  —b^  a/. 


(19) 


The  displacement  in  meters  is  given  by  vaAtb.  Normalizing  this  value  to  the  well 
known  optimum  azimuth  or  Doppler  resolution  A/2,  where  A  =  A/03dB  denotes 
the  real  aperture  of  the  radar,  yields  the  displacement  in  resolution  cells 


vjfh  ~  h )  _  26>3dBfl6(r)A/ 
A/(203dB)  va 


(20) 


Bearing  in  mind  that  the  synthetic  aperture  time  is  much  larger  than  the  pulse  du¬ 
ration  T,  the  displacement  is  seen  to  be  much  more  severe  in  the  Doppler-direction 
(cf.  (11)). 


Fig.  4  illustrates  the  different  displacements  for  range  and  Doppler-direction  as 
well  as  for  typical  airborne  and  spacebome  parameters.  The  pulse  duration  for 
either  case  was  assumed  to  be  T  =  50  /xs.  From  the  solid  curve  in  Fig.  4,  a  target 
becomes  shifted  by  exactly  one  range  cell  for  a  frequency  mismatch  of  A /  =  20 
KHz.  Therefore  the  range  displacement  will  only  be  significant  for  a  rather  large 
discrepancy  of  the  two  local  oscillator  frequencies.  In  contrast,  the  corresponding 
displacement  in  Doppler  direction  will  be  on  the  order  of  10, 000  to  one  million 
pixels  for  this  frequency  offset. 

If,  for  instance,  a  maximum  shift  of  only  one  pixel  is  required,  the  disagreement 
between  the  oscillator  frequencies  must  be  smaller  than  0.1  Hz,  which  is  highly 
unrealistic.  However,  it  is  important  to  note  that  no  de-focussing  in  either  direction 
occurs  due  to  a  constant  frequency  offset,  see  also  [19]. 
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Figure  4  :  Expected  displacement  after  range  and  Doppler  compression  along 
oscillator  frequency  offset,  range  displacement  in  pixels  for  T  — 
50  ps  (solid),  Doppler  displacement  for  Ts  =  0.5  s  (dashed)  and 
Doppler  displacement  for  Ts  =  10  s  (dash-dotted). 


2.2  Linear  oscillator  frequency  drift 

In  addition  to  a  constant  frequency  offset,  in  practise  real  oscillators  drift  slowly 
away  from  their  nominal  center  frequency,  for  instance,  caused  by  environmental 
changes,  such  as  temperature,  pressure  or  magnetic  fields.  Modem  oscillators  com¬ 
monly  compensate  for  such  drifts  by  temperature  compensation.  A  Temperature 
Compensated  Crystal  Oscillator  (TCXO)  typically  contains  a  temperature  compen¬ 
sation  network  to  sense  the  ambient  temperature  and  pull  the  crystal  frequency  to 
prevent  frequency  drift  over  the  temperature  range.  An  Oven  Controlled  Crystal  Os¬ 
cillator  (OCXO)  usually  contains  an  oven  block  (with  temperature  sensor,  heating 
element,  oven  circuitry,  and  insulation)  to  maintain  a  stable  temperature.  OCXOs 
offer  a  great  improvement  in  frequency  stability  compared  to  simple  TCXOs.  Al¬ 
though,  these  techniques  usually  greatly  improve  the  performance,  residual  drifts 
can  still  cause  errors,  particularly  in  radar  applications.  In  space-based  across-track 
interferometry,  for  instance,  very  long  orbital  segments  are  usually  analysed  in  or¬ 
der  to  generate  digital  elevation  models  of  entire  continents.  Even  the  smallest 
clock  drifts  lead  to  measurable  and  compromising  phase  artifacts  [9,  20]. 

The  aim  of  this  section  is  to  investigate  how  frequency  drift  influences  the  per¬ 
formance  of  a  bistatic  SAR,  particularly  since  the  two  oscillators  are  spatially  sepa¬ 
rated.  The  carrier  frequencies  are  assumed  to  drift  linearly  away  from  their  nominal 
value  by 

uji (f)  =  u>i  -f-  ait  and  uj2(t)  =  u>2  -I-  02^1  (21) 
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where  the  a,  denote  the  change  of  rate  of  frequency  or  drift  rates;  the  physical  unit 
is  rad/s.  A  linear  frequency  drift  corresponds  to  a  quadratic  phase  deviation. 

The  quality  of  oscillators  is  commonly  specified  by  a  dimensionless  relative  stabil¬ 
ity  quantity  t  =  A// /,  e.g.,  e  =  1  •  10~8  means  a  frequency  deviation  of  100  Hz 
for  a  10  GHz  X-band  oscillator.  Typical  drift  rates  are  10  to  100  Hz/s  for  TCXOs, 
0.1  to  1  Hz/s  for  higher  quality  OCXOs,  and  smaller  than  0.1  Hz/s  for  the  highest 
quality  oscillators  [21]. 


2.2.1  Fast-time  (range)  processing 

In  order  to  study  the  influence  of  drift  individually,  the  constant  oscillator  frequency 
offset  is  set  to  Acu  =  0,  which  means  that  the  two  oscillators  may  be  perfectly 
synchronized  at  one  point  in  time,  and  then  are  ‘free-running'  afterwards. 


Generalizing  the  received  (and  into  baseband  down-converted)  signal  in  (4)  to  time- 
dependent  oscillator  frequency: 


Srec(f)  =  e-«(t)VIf(,-t"-T)2  e 


,ju>l{t-T)-(t-r) 


rect 


t  —  tn  —  T 


(22) 


inserting  (21)  into  (22)  and  applying  the  filter  function  (5),  leads  to  a  matched  filter 
response  (analogous  to  (6)) 


p(t)  = 


T/2 

-T/2 


■T> 


'gjoi (*-*'— r)2 e-ja2(t-t’)2 e-jn§ t’2  ^ / 


(23) 


where  the  integration  limits  are  restricted  to  one-half  of  the  pulse  duration,  as  in 
(10).  Re-writing  (23)  yields 


p U'j  —  g-JWiTgjTrf  (t-tn-r)2e70i(t-r)2e-ja2t2 
T/2 

f  g— 2j(tt  y  (i— tn— r)+Aai— air^t'gj Aa('2 
-T/2 


(24) 


Since  the  second  order  phase  deviation  is  virtually  always  smaller  than  n/2  over 
the  pulse  duration,  i.e.,  (ai  —  a2)T2  =  A aT2  <§;  7r/2,  the  term  exp(jAaf'2)  can  be 
omitted  in  (24)  and  we  get 

p (£)  _  g-jwiTgjTrf  (t-tn-T)2  ejai(t-r)2  e-ja2t2 

•  sine  f  —  tn  —  r)  +  A  at  —  a\T  j  Tj  .  (25) 
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Assuming  no  range-ambiguities,  the  round  trip  delay  r  is  restricted  to  lie  within 
the  pulse  repetition  interval  r  <  PRI  =  1/PRF  and  the  third  summand  in  the 
argument  of  the  sine-function  can  be  omitted.  The  peak  of  the  sinc-function  occurs 
at  time 


,  nB(tn  +  r) 

0  “  nB  +  AaT' 

which  is  different  from  the  error-free  value  t0  =  tn  +  r  by 


(26) 


A£q  —  t'0  —  to  — 


AaTjtn  +  t)  ^  AaT 
7r  B  +  AaT  tcB 


(27) 


for  large  n.  In  other  words,  the  peak  position  of  the  matched  filter  response  in  the 
range  direction  drifts  very  slowly  away  from  the  correct  location  over  time.  How¬ 
ever,  this  effect  will  only  become  significant  for  unreasonably  long  time  periods. 
For  instance,  it  would  take  more  than  one  hour  to  drift  one  range  cell  (BAt0  =  1) 
when  using  TCXOs,  and  about  100  hours  when  using  OCXOs. 


2.2.2  Slow-time  (Doppler)  processing 


By  omitting  the  drift  (Afo  =  0)  in  (27),  i.e.,  by  sampling  the  signal  (25)  at  slow- 
times  tn  +  r  yields  the  Doppler  history 

SD(tn)  =  e-ia'ire'0lt"e_J'02(tn+T)2  “  e-Mv'A at£g-j2a2<"T,  (28) 

where  a2r 2  is  negligibly  small.  Re-writing  the  time  delay  r  in  (13)  and  (14)  to 


2  Vl  i2\  2 vl.  ,  ,  v\  2  _ 


T  —  -  [  Rb  +  - jj~tb  j  —  —^tbt  +  -j^-t  =  So  +  Sit  +  s2t  , 


(29) 


cRb  cRb 

the  Doppler  signal  in  (28)  becomes 

SD(f)  =  g-j(u’lal+2°2«o)te-i(‘*'iS2-l-2a2Si-Aa)t2e-j2a2S2t3  (30) 

where  for  notational  simplicity  tn  was  replaced  by  t  and  the  constant  term  omitted. 
After  application  of  the  reference  function  sIe((t)  =  exp(ju>is2t2)  rect (t/Ts),  the 
point  spread  function,  cf.  ( 1 8),  yields 

T./2 

pD(t)  =  J  e-j(u,isl+2“2so)(t-t,)e-j(wi»2+2a2Si— Aa)(t— «,2)e— i2o2S2(t— t')3gj“l«2 t'2^' 


-T,/  2 
-j(ui 
T./ 2 


_  g— j(wisi+2a2So)ig— j(‘*’iS2+2a2Si— Aa)(2g— j2a2S2t3 


/ 


gj(uisi+2a2So+2(a;iS2+2a2Si  -  Aa)t+6a.2S2t2)t' 


-T./2 

.  gj(2a2si  -Aa+6a2S20t'2g-J2a2«2«'3  ^  Q  j  j 

Two  conditions  must  be  fulfilled  to  obtain  a  sharply  focussed  point  spread  function: 
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Table  1:  Second  order  term  Ts2Aa  for  oscillator  drifts 


TCXO 

Aa  =  10  Hz/s 

ocxo 

Aa  =  0.1  Hz/s 

oocxo 

Aa  =  0.01  Hz/s 

Airborne:  Ts  =  10  s 

1000 

10 

1 

Space-borne:  Ts  =  Is 

10 

0.1 

0.01 

1.  |2a2s2|  Tg  <  7r/2  =►  M  <  *£*  = 

2.  |  Aa  —  2a2si  —  6a2s2f|  Ts2  <  7r/2  =*>•  |Aa|  < 

While  the  first  inequality  is  always  fulfilled  because  si,s2  =  0,  the  second  one 
depends  on  the  maximum  synthetic  aperture  time  Ts.  From  Table  1  it  is  evident  that 
even  small  oscillator  frequency  drifts  may  cause  a  degradation  and  widening  of  the 
point  spread  function  in  the  Doppler  direction.  This  loss  of  Doppler  resolution  will 
be  significant,  particularly  for  airborne  systems  with  larger  synthetic  aperture  times 
compared  to  their  space-borne  counterparts.  It  is  clear  that  special  SAR  modes  with 
higher  dwell  times,  such  as  spotlight  mode,  will  increase  this  problem  to  the  point 
that  the  focussing  behavior  may  be  completely  lost. 

2.3  Random  oscillator  phase  noise 

A  further  error  source  for  bistatic  imaging  is  the  unavoidable  physical  phenomenon 
of  phase  noise  or  jitter  of  the  oscillators.  Jitter  describes  the  frequency  instabilities 
due  to  random  noise  processes  in  any  real  electronic  components  [22],  The  effect 
of  multiplicative  phase  noise  on  bistatic  radar  performance  has  been,  for  instance, 
investigated  in  [23]. 

In  contrast  to  deterministic  oscillator  (low-frequency)  noise  types,  such  as  fre¬ 
quency  offset  or  drift,  which  result  in  geometric  distortions  (e.g.,  loss  of  resolution 
and  increase  of  peak  sidelobes),  wideband  phase  errors  spread  out  energy  widely 
across  the  impulse  response  or  point  spread  function.  Rather  than  causing  a  smear¬ 
ing  of  the  mainlobe  or  increasing  of  the  first  sidelobes.  the  phase  jitter  contributes  to 
an  increased  sidelobe  level  over  the  entire  sidelobc  range,  which  results  in  contrast 
degradation  and  noise  in  the  image.  A  measure  of  such  an  error  is  the  integrated 
sidelobe  ratio  (ISLR).  The  ISLR  is  the  ratio  of  all  energy  in  the  sidelobes  to  the 
energy  in  the  mainlobe  of  the  impulse  response.  In  its  common  definition,  the  di¬ 
vision  between  mainlobe  and  sidelobe  is  the  location  of  the  first  zero  crossing  of 
the  noise-free  impulse  response.  The  significance  of  the  integrated  sidelobe  en¬ 
ergy  depends  not  only  on  the  nature  of  the  sidelobes,  but  equally  important  on  the 
distribution  of  scatterers  throughout  the  scene.  If  the  scene  contains  mainly  a  few 
sparsely  distributed  bright  scatterers,  then  the  ISLR  is  less  important  than  the  peak 
sidelobe  ratio  (PSLR),  which  describes  the  amount  of  energy  in  the  highest  sidelobe 
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compared  to  the  mainlobc.  High  peak  sidclobcs  appear  as  false  targets  and  hinder 
detection  of  closely  located  weaker  targets.  In  contrast,  a  low  ISLR  is  important  for 
homogeneous  scenes  with  evenly  distributed  weaker  scatterers,  where  a  high  ISLR 
will  degrade  contrast  by  filling  in  low  returns  and  shadow  areas  [24], 

2.3.1  Wideband  (white)  phase  noise 

In  order  to  quantify  the  resulting  ISLR  for  bistatic  SAR,  we  start  with  the  com¬ 
pressed  pulse  (23),  where  the  drift  terms  are  replaced  by  the  corresponding  phase 
noise  components  (t)  and  £2(t)  (cf.  also  (2)) 

T/2 

P(t)  —  e-Ja'lT^;  /  <j*T  (t-t'-tn-r)^E1(t-t'-T)e-jE2(t-t')e-jn§t^  dt,  (32) 
-T/2 

The  capital  letters  in  (32)  are  used  to  distinguish  stochastic  processes  from  their  cor¬ 
responding  realisations.  The  phase  noise  terms  are  modeled  as  stationary  stochastic 
Gaussian  processes  E (t)  with  expectation  zero  and  auto-covariance  function 

OO 

i(t)  =  J  c3i  Ei(u)  exp  (jut)  du  with  ce,h,(0)  =  of,  (33) 

—  OC 

where  Cs^iu)  denotes  the  spectral  density  function  or  the  power  spectrum.  The 
spectral  density  function  is  an  important  performance  measure  which  is  usually  sup¬ 
plied  by  the  manufacturer  of  oscillators.  For  a  conventional  monostatic  SAR  where 
the  same  oscillator  is  used  for  up-  and  down-conversion,  32  equals  Ei  (cf.  section 
2.3.3),  whereas  in  the  bistatic  case,  two  independent  oscillators  are  used.  Hence, 
E2  and  Ej  are  statistically  independent.  More  details  and  background  information 
about  the  statistical  frequency  characteristics  of  precision  clocks  and  oscillators 
can,  for  instance,  be  found  in  [22,  25,  26]. 

As  a  function  of  two  random  processes,  the  impulse  response  function  (32)  is  again 
a  stochastic  process.  For  statistically  white  phase  noise,  its  statistical  properties 
have  been  determined  in  terms  of  the  expectation  function 

E  P(t)  =  e  L2-2' sine  (n B(t  —  tn  —  r))  (34) 

and  the  second  moment  function 
( 1  - 

E  \P(t)\2  =  y/2 - - — - —  +  e~^+a^  sine2  (n B(t  —  tn  —  r))  (35) 
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Table  2:  Degradation  of  ISLR  due  to  phase  noise. 


IISLR  [dB] 

cr  [rad] 

-5 

0.56 

-10 

0.32 

-15 

0.18 

-20 

0.10 

-25 

0.056 

-30 

0.032 

in  appendix  A.  Therefore,  the  variance  of  P(t)  can  be  computed  as 

o2p{t)  =  E \P(t)\2  -  |E  P(t) |2  =  ^  (l  -  e-^+^)  ,  (36) 

which  denotes  the  power  of  the  random  part  of  the  stochastic  impulse  response,  and 
therewith  can  be  used  to  determine  the  sidelobe  power  levels  caused  by  the  phase 
jitter. 

Integrating  (36)  over  the  entire  pulse  width  T,  and  using  exp(— cr2)  «  1  —  cr2  for 
small  phase  noise  terms,  it  becomes  evident  that  the  power  level  is  raised  over  the 
entire  compressed  pulse  length  by  the  sum  of  the  two  phase  noise  variances.  Hence, 
the  increase  (i.e.,  degradation)  of  the  ISLR  (IISLR)  compared  to  the  undisturbed  or 
phase  noise-free  case  is 

OO  OO 

IISLR  =  V2( a2  +  a2)  =  V2  j  C2lSl  (/)  +  CE2h2(/)  df  =  4^2  J  C==(/)  df, 

— oo  0 

(37) 

where  it  was  assumed  that  identical  oscillator  types  with  o\  —  a2  are  used.  The 
above  analysis  is  analogously  applicable  to  the  Doppler  direction,  where  the  pulse 
parameters  in  (32)  are  replaced  by  the  Doppler  rate  in  (17),  and  the  pulse  duration 
T  is  replaced  by  the  synthetic  aperture  time  Ts. 

Table  2  lists  the  allowable  phase  noise,  more  specifically  the  constraint  on  the  stan¬ 
dard  deviation  or  root  means  square  (rms)  value  a,  in  order  not  to  exceed  a  certain 
degradation  in  ISLR  compared  to  the  noise-free  case.  A  maximum  ISLR  degrada¬ 
tion  due  to  phase  noise  of  IISLR  =  —  20  dB  is  a  typical  requirement  for  imaging 
applications  [24], 

To  illustrate  the  impact  of  phase  noise  graphically,  the  red  curve  in  Fig.  5  shows 
the  effect  of  a  a  =  0.16  rms  phase  jitter  (IISLR  =  -16  dB)  on  the  impulse  response 
function  in  a  bistatic  scenario.  In  contrast,  the  case  of  cr  =  0.04  is  plotted  in  blue 
(IISLR  =  -28  dB),  which  is  virtually  indistinguishable  from  the  error-free  case. 
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Figure  5  :  Impulse  response  function  in  the  presence  of  white  phase 
noise  with  o  =  0.04  (blue)  and  a  =  0.16  (red). 


2.3.2  Realistic  (colored)  oscillator  phase  noise 

In  practical  circumstances,  the  aforementioned  low-frequency  and  high-frequency 
components  of  the  phase  noise  will  always  occur  simultaneously.  In  order  to  predict 
the  performance  degradation  due  to  phase  noise,  e.g.,  the  maximum  achievable 
coherent  processing  time,  the  single  sideband  (SSB)  power  density  spectrum  £(/) 
measured  around  the  particular  oscillator  output  frequency  /0  can  be  used.  This 
figure  is  commonly  provided  by  the  device  manufacturers.  The  SSB  and  the  spectral 
density  in  (33)  are  related  by  the  multiplication  factor  M,  which  is  needed  to  up- 
convert  /0  to  the  transmit  frequency  fc: 


Os(/)  =  2M2C(f)  =  2  (^j  £(/).  (38) 

As  examples.  Table  3  shows  the  provided  SSBs  for  three  commercially  available 
oscillators,  the  58503  from  Symmetricom  Inc.,  the  high  quality  oscillator  405681 
from  Wenzel  Inc.,  and  the  ultra  low  phase  noise  oscillator  from  Oscilloquarz  Inc. 
Inserting  these  values  into  the  white  noise  formula  (38)  and  integrating  the  result 
over  the  entire  stated  frequency  range  according  to  (37)  yields  rough  estimates  for 
the  phase  noise  rms  expected  from  these  three  oscillators.  The  calculated  variance 
for  the  highest  quality  OSC  is,  for  instance,  about  37r  and  for  the  others  even  larger. 
Evidently,  these  values  are  much  greater  than  one  radian.  Since  the  white  noise 
assumption  is  not  valid,  the  small  noise  approximation  in  (37)  cannot  be  applied, 
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Table  3:  Measured  SSBs  for  three  commercial  oscillators. 


Type 

Symmetricom  58503B 

Wenzel  405681 

Oscilloquarz  (OSC) 

/o  [MHz] 

10 

10 

5 

/[Hz] 

A/)  [dB] 

0.001 

5 

-15 

-40 

0.01 

-25 

-45 

-70 

0.1 

-55 

-75 

-100 

1 

-85 

-105 

-130 

10 

-125 

-135 

-145 

100 

-135 

-160 

-152 

1000 

-140 

-176 

-155 

10000 

-145 

-176 

-157 

i.e..  the  I1SLR  cannot  be  specified  by  this  formula.  Nevertheless,  is  it  obvious  that 
all  oscillator  types  exceed  substantially  the  requirement. 

In  addition,  these  estimated  values  represent  the  phase  jitter  variance  or  standard 
deviation  based  on  a  very  large  time  base.  In  the  considered  case,  it  is  the  phase 
rms  over  T  =  1/ fmin  =  1000  s.  In  practice,  however,  the  time  interval  over  which 
any  decorrelation,  i.e.,  random  phase  variation,  will  effect  the  impulse  response 
function  is  much  smaller.  For  bistatic  SAR,  this  time  is  identical  to  the  synthetic 
aperture  time  Ts,  which  typically  ranges  from  about  one  second  in  space-based 
SAR,  to  a  few  seconds  in  airborne  stripmap  SAR,  and  up  to  hundreds  of  seconds 
in  wide-range  spotlight  mode  SAR  applications  [27].  For  finite  Ts,  the  rms  phase 
value  er  has  been  calculated  in  appendix  B;  it  can  be  determined  via  the  integral 


Reducing  the  lower  integration  bound  to  about  1  / (2 Ts)  was  found  to  give  more  ac¬ 
curate  results  than  the  limit  l/Ts  stated  in  some  publications,  e.g.,  [23],  The  black 
curve  in  Fig.  6  shows  the  IISLR  (10log(cr2))  in  (39)  for  the  OSC  oscillator  over 
increasing  synthetic  aperture  time  Ts,  which  was  varied  from  0.1  s  to  50  s.  The 
point  of  intersection  with  the  -20  dB  limit  for  the  IISLR  indicates  a  maximum  Ts  of 
about  1 0  seconds. 

However,  even  the  rms  phase  in  (39)  does  not  reveal  the  full  impact  of  jitter  on 
the  point  spread  function.  This  is  illustrated  in  Fig.  7,  where  the  undisturbed  point 
spread  function  is  shown  along  with  several  realisations  that  include  phase  noise, 
which  are  simulated  via  OSC’s  spectral  density  in  Table  3.  Although  Fig.  6  would 
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indicate  a  relatively  large  IISLR  of  -10  dB  for  the  chosen  time  Ts  =  30  s,  almost 
no  random  raise  in  the  sidelobes  further  away  from  the  mainlobe  are  recognisable 
in  Fig.  7  (a).  In  contrast.  Fig.  5  shows  a  significant  increase  in  the  sidelobes  for  an 
evidently  smaller  IISLR  of  -16  dB. 


Figure  6  :  IISLR  of  OSC’s  oscillator  for  increasing  synthetic  aperture  time, 
upper  curve  based  on  eq.  (39),  lower  curve  based  on  eq.  (42). 

This  seeming  contradiction  can  be  explained  with  the  fact  that  (39)  does  not  dis¬ 
tinguish  between  the  different  effects  of  low-frequency  and  high-frequency  com¬ 
ponents  of  the  phase  noise.  Low-frequency  phase  noise  components  will  lead  to 
geometric  distortions  (cf.  Section  2.2)  but  no  increase  in  sidelobe  power,  which 
in  fact  is  exclusively  caused  by  the  high-frequency  components,  see  Section  2.3. 
Fig.  7  (b)  confirms  the  randomly  increased  peak  sidelobes  due  to  the  low-frequency 
components. 

Therefore,  a  much  more  meaningful  metric  to  describe  the  overall  phase  noise  im¬ 
pact  must  remove  the  low-frequency  (slowly-time  varying)  components  before  de¬ 
termining  the  IISLR.  For  instance,  we  can  model  these  low-frequency  components 
as  a  low-order  polynomial  function  in  time.  In  most  cases,  an  order  of  three  for  the 
phase  function  (corresponding  to  quadratic  frequency  terms)  is  sufficient: 

—(f)  =  bo  +  b\t  +  £>2^  *”  +  ,J  +  Y(f),  (40) 

where  T(f)  represents  the  remaining  high-frequency  phase  noise.  A  similar  ap¬ 
proach  has  been  used  to  predict  the  baseline  accuracy  between  platforms  in  bistatic 
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(a)  Sidelobe  region  around  -50  dB.  (b)  Mainlobe  region. 

Figure  7:  Ideal  point  spread  function  (black)  along  with  ten  realisations  including 
phase  noise  with  spectral  density  in  the  right-hand  column  of  table  3. 


interferometry  [28],  in  which  a  direct  link  between  the  two  radars  with  a  maximum 
decorrelation  time  in  the  order  of  the  PRI  was  proposed.  On  the  other  hand,  we  are 
interested  in  indirect  synchronisation  over  significantly  longer  integration  times. 

As  mentioned  before,  estimates  of  the  coefficients  b  =  [60,  -  -  • ,  b3]'  will  give  insight 
into  the  geometric  distortions  of  the  point  spread  function,  and  the  variance  of  T 
will  describe  the  radiometric  distortions,  i.e.,  the  increase  in  ISLR.  Taking  both  ef¬ 
fects  into  account  simultaneously  will  then  yield  the  maximum  allowable  coherent 
processing  or  synthetic  aperture  time.  Since  the  phase  noise  process  <!>(£)  is  as¬ 
sumed  to  be  normally  distributed,  all  estimated  coefficients  6*  are  shown  (annex  C) 
to  be  Gaussian  random  variables  with  zero  mean  and  variance 


oo 

<=4  J Chaif)  \Vi(f)\2df  i  =  0 . 3, 

0 

(41) 

where  V*(/)  is  the  Fourier-transform  of  Vi(t)  defined  in  (C.I3). 

Similarly,  the  variance  of  the  remaining  high-frequency  phase  noise  components 
are  calculated  in  annex  C  as 

OO 

4=4/ CehM)  \Mf)\2df, 

0 

(42) 

where  it  can  be  shown  that  V*(/)  and  A(f)  have  high-pass  character  with  increasing 
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cut-off  frequencies  as  larger  polynomial  orders  in  (40)  are  chosen.  Therefore,  the 
larger  the  order  the  lower  the  integrated  power  under  the  spectral  density  function 
(cf.  (39))  and  the  smaller  the  variance  of  the  remaining  high-frequency  phase  noise. 

In  summary,  the  derived  statistics  of  the  polynomial  coefficients  and  the  variance 
of  the  remaining  high-frequency  phase  noise  can  be  used  to  predict  the  geometric 
and  radiometric  distortions  of  the  point  spread  function.  The  constant  parameter  b0 
is  unimportant  and  the  linear  term  bx  is  circumstantial  because  it  introduces  only 
a  shift  of  the  function,  which  can  also  be  seen  in  Fig.  7.  Reduced  resolution  and 
increased  peak  sidelobes  are  due  to  b2  and  63  only.  As  a  rule  of  thumb,  these  pa¬ 
rameters  must  not  exceed  the  following  thresholds  (compare  (3 1 )) 

N  *  ^  and  N<^-  <43) 

The  magnitude  of  the  normally  distributed  estimates  b2  and  b:i  are  governed  by  the 
probability  density  function  (pdf)  f^ih)  =  2/^.(6j)u(6j)  where  u(bi)  denotes  the 
unity  step  function  [29,  30].  The  mean  value  of  each  estimator  is  give  as 

E  ^  =  yf  a'bi  f°r  *  =  °’  ‘  ‘  ‘  ’ 3’  (44) 

and  hence  depends  only  on  the  variance  in  (41).  Figure  8  compares  these  mean 
values  of  the  polynomial  coefficients  with  the  thresholds  in  (43)  for  increasing  syn¬ 
thetic  aperture  times  Ts.  The  two  curves  intersect  with  the  thresholds  at  roughly 
11  s  for  63  and  32  s  for  62,  respectively,  where  the  smaller  value  obviously  sets 
the  upper  bound  for  the  maximum  Ts.  For  confirmation  of  the  theoretical  results, 
the  mean  value  of  ten  simulated  coefficients  for  each  T,  are  superimposed  showing 
almost  perfect  agreement  with  the  theoretical  results. 


The  corresponding  remaining  IISLR  after  subtraction  of  these  estimates  is  plotted 
in  Fig.  6  (red  curve),  where  it  it  is  evident  that  the  sidelobc  levels  are  now  always 
below  the  requirement  of  -20  dB  for  T,  <  50  s. 


Figure  8  does  not  shed  light  on  how  often  we  expect  that  geometric  distortions  will 
occur.  In  fact,  by  comparing  only  the  mean  value,  we  know  that  about  50  %  of  the 
time  the  values  will  be  larger  than  the  threshold,  and  the  other  50  %  will  be  smaller. 
Therefore,  we  may  calculate  the  probability  that  |6j|  is  larger  then  the  corresponding 
threshold, 


P 


l-p  {n< 

2  [  1  -  P  i  dk 


(45) 
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where  $(•)  is  the  cumulative  distribution  function  of  the  standard  normally  dis¬ 
tributed  random  variable  dj  (error-function),  i.e.,  with  zero  mean  and  unity  variance. 
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Figure  8  :  Mean  values  of  the  polynomial  coefficients  with  corresponding 
thresholds  versus  synthetic  aperture  time. 
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Figure  9  shows  the  corresponding  probabilities  for  the  quadratic  and  cubic  coeffi¬ 
cients.  If,  for  instance,  one  is  only  willing  to  accept  some  geometric  distortion  one 
out  of  ten  times  (i.e.,  10  %),  e.g.,  one  out  of  ten  flight  lines  during  an  experiment, 
then  the  coherent  processing  time  must  be  limited  to  about  Ts  =  5  s. 

It  was  found  by  simulation  that  these  conditions  are  rather  conservative  and  that  the 
thresholds  can  in  many  case  be  relaxed  to 

N  <  ^  and  \b3\  <  (46) 

for  which  the  probabilities  are  also  plotted  in  Fig.  9.  For  the  same  number  of  al¬ 
lowed  oversteppings,  i.e.,  10  %,  the  maximum  allowed  synthetic  aperture  time  is 
now  about  15  s. 

Although  Fig.  9  gives  an  indication  of  how  often  the  thresholds  are  exceeded,  it 
does  not  show  by  how  much,  i.e.,  the  actual  impact  on  the  point  spread  function. 
In  other  words,  if  the  thresholds  are  only  exceeded  by  a  relatively  small  amount, 
the  geometric  distortions  (such  as  wider  mainlobe  or  non-symmetric  sidclobes)  are 
expected  to  be  negligible.  An  indicator  of  such  behavior  can  be  seen  in  Fig.  8,  where 
the  mean  values  of  b2  and  b:i  arc  shown  to  stay  quite  close  to  the  thresholds  for 
increasing  Ts  behind  the  intersection  point.  This  observation  is  also  confirmed  by 
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the  simulation  in  Fig.  7  (b),  where  only  minor  distortions  (1  to  2  dB  increase  of  peak 
sidelobe)  can  be  recognized  even  for  the  relatively  long  Ts  of  30  s.  In  particular, 
the  point  spread  function  appears  perfectly  symmetric  although  a  significant  cubic 
phase  term  was  predicted  in  Fig.  8  (two  bottom  curves).  A  further  reduction  of 
these  sidelobe  contributions  (distortions)  due  to  the  low-frequency  phase  noise  may 
be  achieved  by  using  weighting  functions  (e.g.,  Hamming)  during  the  convolution. 


Figure  9  :  Probability  that  the  polynomial  coefficients  are  larger  than  the  re¬ 
quired  threshold. 


As  a  consequence,  a  performance  analysis  based  on  simulations  may  be  preferable 
in  some  cases  compared  to  pure  theoretical  predictions  via  those  metrics  derived 
above.  Simulations  can,  for  instance,  be  used  to  find  those  values  of  Ts  for  which 
the  3  dB-width  of  the  mainbeam  increases  by  a  given  amount,  or  for  which  the 
peak  sidelobe  has  been  doubled,  etc.  Such  simulations  can  be  performed  by  simply 
using  the  given  (measured  or  provided)  spectral  density  function  of  any  arbitrary 
oscillator  as  the  sole  input. 

Beyond  the  possibility  of  predicting  and  avoiding  the  performance  degradation  of 
bistatic  SAR  resolution  in  the  presence  of  particularly  severe  phase  noise  by  simply 
using  higher  quality  oscillators,  there  exists  the  opportunity  to  estimate  and  remove 
the  low-frequency  (slowly  varying)  phase  error  prior  to  imaging.  The  applicable 
techniques  are  identical  to  those  used  for  ‘auto-focussing’,  i.e.,  motion  compensa¬ 
tion  based  on  the  measured  data  alone  [31,  32].  The  suitability  of  some  of  these 
techniques  will  be  explored  in  phase  II  of  this  project. 
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Another  important  problem  not  addressed  so  far  is  the  impact  of  motion  and/or  vi¬ 
bration  of  the  oscillator  to  its  phase  noise  characteristics.  It  is  well  demonstrated 
that  oscillator’s  nominal  frequency  reacts  sensitively  to  motion  and  accelerations, 
which  is  usually  provided  as  the  parameter  ‘g-sensitivity’  by  the  manufacturer.  Os¬ 
cillators  with  small  /0  suffer  more  from  this  deficiency  because  of  their  increased 
physical  size  compared  to  higher  frequency  oscillators.  In  contrast,  lower  frequency 
oscillators  commonly  possess  a  much  better  (lower)  phase  noise  characteristic. 
Some  preliminary  measurements  are  introduced  in  section  6.3.  However,  exten¬ 
sive  laboratory  experiments  are  being  executed  during  the  writing  of  this  report  and 
will  also  be  addressed  in  phase  II. 

2.3.3  Comparison  to  monostatic  case 

It  is  important  to  realize  that  the  requirements  on  oscillator  quality  with  respect  to 
phase  noise  are  much  more  stringent  in  bistatic  than  in  monostatic  SAR  due  to  the 
use  of  two  separate  oscillators.  In  monostatic  SAR  the  same  oscillator  is  used  to 
up-  and  down-convert  the  signal  for  each  pulse.  Although  the  oscillator  frequency 
varies  along  azimuth-time  (according  to  its  power  spectrum),  the  received  signals 
can  maximally  decorrelate  over  the  round-trip  travel  time  (PRI)  of  each  pulse.  This 
can  be  seen  by  replacing  the  combined  spectral  density  function  CEE(u;)  in  (37) 
with  its  monostatic  counterpart.  Consider  the  phase  term  difference  in  (32)  for  the 
monostatic  case  E(£)  =  Ex(t  —  r)  —  Ej  (i)  for  which  the  covariance  function  can  be 
calculated  as 

cHH(0  =  EE  (t  +  t')E*(t) 

=  E  (Ex(f  +  t'  -  t)  -  E1(t  + 1'))  (Sj(t  -  r)  -  E,(f))* 

=  2 C£l3l  it')  -  c=lSl  (f'  +  r)  -  cElEl  ( t '  -  t).  (47) 

Using  the  shift  property  of  the  Fourier-transform  yields  for  the  spectral  density 
function 

CeeM  =  2ChiSi  M  -  CejEi  Me*"  -  CSlSl  (u)e~^ 

=  2(1  -  cos(urr))CElHlM  =  4sin2  (y)  (48) 


Equation  (47)  reveals  a  kind  of  high-pass  filtering  to  the  spectral  density  function 
because  low-frequency  components  are  suppressed  with  a  sin2-charactcristic  below 
frequencies  of  1/r.  Integration  of  (48)  over  this  frequency  range  yields 
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which  corresponds  to  an  monostatic  IISLR  of  —28  dB  for  the  Symmetricom,  —40 
dB  for  the  low  noise  oscillator  from  Wenzel,  and  about  —52  dB  for  the  ultra-low 
phase  noise  oscillator  form  OSC. 
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3  Bistatic  observation  time 


One  parameter  of  major  importance  for  bistatic  SAR  is  the  time  period  for  which  a 
scatterer  on  a  particular  location  on  the  ground  is  observed  simultaneously  by  both 
radars.  In  order  to  calculate  this  time,  the  individual  time-dependent  range  vec¬ 
tors  (1)  for  the  transmitter  and  receiver  are  transformed  into  the  antenna  coordinate 
system  by 


for  i  =  1,2,  where  denotes  the  constant  course  angle  defined  in  appendix  E.3. 
<1 \{t)  and  ty(t)  are  the  azimuth  and  elevation  (depression)  angles  (appendix  E.4) 
which  are  in  general  adjusted  over  time,  depending  on  the  particular  SAR  mode 
used: 

•  Stripmap  or  squinted  stripmap  mode:  <&i(t)  =  and  ^i(t)  =  where 

o  and  tfi|0  are  the  squint  angle  and  depression  angle  to  the  center  of  the  scene 
at  time  t  =  0. 

•  Spotmode:  $,:(£)  and  'l'j(t)  are  updated  from  pulse  to  pulse  so  that  the  antenna 
is  always  pointing  at  center  of  the  scene  (r  =  0). 

It  is  obvious  that  the  bistatic  observation  time  can  almost  be  arbitrarily  chosen  in 
spotmode  operations,  as  long  as  the  scatterer  is  within  the  overlap  region  of  both  an¬ 
tenna  footprints.  For  stripmap  operations,  however,  where  the  beam  pointing  of  the 
antenna  is  kept  fixed  over  the  measurement  interval,  the  bistatic  observation  time  is 
principally  limited.  Depending  on  the  flight  trajectories  of  both  platforms,  the  max¬ 
imum  observation  time  may  vary  significantly  for  different  areas  on  the  ground.  In 
order  to  assess  whether  a  particular  point  on  the  ground  is  simultaneously  ‘seen’  at 
a  certain  instant  in  time,  the  directional  cosine  u,  and  sine  of  both  radars  must  be 
computed: 


Ui(t,r) 


Vi(t,r)  = 


(51) 


where  ux  =  [1,0,0]'  and  uz  =  [0,0, 1]'.  The  bistatic  observation  time  7"b  is  then 
given  as  the  integral 


OO 


(52) 


—  OO 
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where  the  indicator  functions  are  defined  as 


{1  if  \ui(t,r)  -  uoti|  <  U3dB,i/2 

\vi(t,r)  -  vto.il  <  V3dB,i/2  (53) 

0  elsewhere. 


u0,i  and  v0j  are  the  directional  cosine  and  sine  corresponding  to  the  azimuth  and  el¬ 
evation  angles  and  ty0,;  of  both  platforms  to  the  scene  center  at  time  origin.  The 

antenna  beamwidth  in  either  direction  are  herein  denoted  by  the  subscript  ‘3dB\ 

By  inserting  (1)  into  (51),  it  is  possible  to  solve  the  inequality  of  (53)  analytically. 
The  derivation  is  exemplified  below  for  the  directional  cosine  where,  for  conve¬ 
nience,  the  transmitter/receiver  subscript  i  is  omitted: 


u(t,r) 


y^LJJlo  -  r)  +  y^L^vt 

\/ILEo  -  ill2  +  2CSo  -  r)'vt  +  ||v||2*2 


<  »7> 


(54) 


where  t?  =  u0  +  u^/2,  which  can  be  re-written  as 


(do  +  di  t)2 
C2  t2  +  C\  t  +  Cq 


<  rj2  (bo  +  b\  t  +  62  t 2) 
=  0, 


with 


(55) 


a0  =  y^LasiEo  -  r) 

bo  =  ILSo  “ill2 
bi  =  2 (Rq  -  r)'v 

b2  =  INI2,  (56) 

and  Co  =  d^  —  r]2b0,  C\  =  2d0di  —  rfbx  and  C2  =  df  —  r/2^-  The  solution  of  (55)  is 
given  as 

*1,2  <  -p-  T  \  -—  +  •%  for  c2  ±  0,  (57) 

2c2  y  c2  4c^ 

Depending  on  the  sign  of  c2,  different  cases  for  the  solution  must  be  distinguished. 
According  to  Fig.  10,  there  is  either  no  real  solution  if  the  parabola  does  not  cross 
the  x-axis,  or  two  real  solutions  in  the  specified  gray  regions. 
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Figure  10  :  Regions  of  possible  solutions  for  the  bistatic  observation  time. 


Comparing  the  minima  and  maxima  of  the  individual  solutions  for  u  and  v  for  both 
transmitter  and  receiver  yields  the  maximum  observation  time  for  each  scatterer  on 
the  ground: 

Tb(r)  =  max  (0,  tmax(r )  -  tmin(r)) ,  (58) 

where 

tm-m  =  max^i1,*!2,^1,^2) 

fraax  =  min  (59) 

and  where  ‘max(:c,  y,  z )’  means  the  largest  of  each  value  x,  y  and  2. 

Figure  1 1  demonstrates  the  variation  of  bistatic  observation  time  Tb(r)  given  by 
(58)  for  each  pixel  on  the  ground  for  different  bistatic  angles  between  transmitter 
and  receiver.  In  this  case,  the  transmitter  is  chosen  to  be  an  aircraft  with  parameters 
listed  in  Table  4.  The  receiver  is  located  on  top  of  a  stationary  tower  with  parame- 


Table  4:  Airborne  transmitter  parameters. 


Altitude 

III  =  5000  m 

Velocity 

Uj  =  132  y^.  m/s 

Elevation  angle  to  scene 
centre  at  t  =  0 

$0,1  =  10° 

Squint  angle  to  scene 
centre  at  t  =  0 

O 

O 

II 

o' 

Antenna  pattern 

sine 

Azimuth  beamwidth 

^3dB,l  =  4° 

Elevation  beamwidth 

^3dB,l  =  2.5° 

ters  listed  in  Table  5,  where  the  receiver  squint  (azimuth)  angle  ^0,2  is  varied,  see 
appendix  F.  ^0,2  =  0°  means  the  aircraft  is  flying  behind  the  tower  and  both  an¬ 
tennas  are  pointing  in  the  same  direction  (i.e.,  a  quasi-monostatic  geometry).  In 
contrast,  ^0,2  =  0°  indicates  that  the  aircraft  is  flying  in  front  of  the  radar  perpen¬ 
dicular  to  the  receiver  antenna  pointing  direction,  resulting  in  a  forward-scattering 
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Table  5:  Tower  receiver  parameters. 


Altitude 

H2  =  20  m 

Velocity 

v2  =  0  m/s 

Elevation  angle  to  scene 
centre  at  t  =  0 

$0,2  —  0.5° 

Antenna  pattern 

sine 

Azimuth  beamwidth 

U3dB,2  =  45° 

Elevation  beamwidth 

^3dB,2  =  10° 

geometry.  It  is  evident  that  the  maximum  observation  time  is  relatively  large  ( Ts  « 
12  to  16  s)  because  of  the  shallow  elevation  angle  of  the  transmitter,  which  results 
in  large  distances  to  the  scene  on  the  ground.  A  larger  distance  means  an  increased 
synthetic  aperture  time  for  a  given  azimuth  beamwidth.  Also  recognizable  are  the 
different  regions  on  the  ground  which  are  bistatically  observable  due  to  the  change 
in  receiver  location  and  orientation.  Of  course,  much  more  complex  illumination 
shapes  are  usually  possible  when  the  receiver  is  also  moving. 
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y-axis  (m)  y-axis  (m)  y-axis  [m| 


(a)  Quasi-monostatic  geometry,  bistatic  an-  (b)  Bistatic  angle  =  45°. 

gle  =  0°. 


(c)  Bistatic  angle  =  90°. 


(d)  Bistatic  angle  =  135°. 


(e)  Bistatic  angle  =  150°.  (f)  Forward-looking  geometry,  bistatic  angle 

=  180°. 

Figure  11:  Bistatic  observation  time  of  ground  range  scatterer.  Stationary  re¬ 
ceiver  and  airborne  transmitter  for  varying  bistatic  angles  ^0,2* 
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4  Bistatic  SAR  performance  analysis 


The  aim  of  this  section  is  to  derive  tractable  mathematical  tools  to  determine  the 
bistatic  SAR  resolution  for  individual  points  on  the  ground  for  arbitrary  transmitter 
and  receiver  flight  geometries.  The  term  resolution  in  this  context  means  the  ability 
to  separate  close  targets  in  different  areas  on  the  ground.  A  similar  approach  was 
applied  for  a  parasitic  SAR  system  in  [15]. 

These  tools  will  allow  performance  comparisons  of  different  flight  geometries  with¬ 
out  the  necessity  of  SAR  processing,  i.e.,  the  complex  and  sometimes  time-consuming 
computation  of  the  point  spread  function,  e.g.,  [17].  They  could  also  be  used  in  an 
opposite  way  during  the  flight  planning  phase  of  operations  to  determine  the  neces¬ 
sary  flight  geometry  (and  velocities)  for  a  required  minimum  ground  resolution  in 
a  particular  region  of  interest. 

It  is  important  to  note  that  the  analysis  herein  is  based  on  a  fully  three-dimensional 
scenario,  i.e.,  it  is  not  restricted  to  two-dimensional  geometries  as  is  sometimes 
done  in  the  literature.  Further,  constant  velocities  for  the  radar  platforms  are  only 
chosen  for  clarity  of  the  derivation.  Accelerations  or  non-constant  altitudes  can 
easily  be  included  in  the  analysis. 


4.1  Preliminaries:  Doppler-bandwidth, 
Doppler-gradient  and  the  k-space 

In  order  to  determine  SAR  resolution,  it  is  crucial  to  understand  the  physical  phe¬ 
nomena  describing  the  composition  of  the  synthetic  aperture  and  therewith  allow¬ 
ing  the  imaging  of  the  desired  reflectivity  with  high  resolution.  This  is  even  more 
important  as  the  general  bistatic  case  is  usually  much  more  complex  than  the  rela¬ 
tively  simple  monostatic  case.  For  the  monostatic  case,  for  instance,  it  is  common 
in  the  literature  to  use  the  Doppler-bandwidth  (the  maximum  spread  of  Doppler 
frequency)  of  a  scatterer  on  the  ground  to  calculate  the  corresponding  azimuth 
(cross-range  or  Doppler)  resolution.  In  this  subsection  we  want  to  show  that  the 
idea  of  bandwidth  can  only  be  applied  to  the  monostatic  non-squinted  case.  In¬ 
stead  of  Doppler-bandwidth,  the  shape  of  the  point  spread  function  (and  therewith 
the  resolution)  is  determined  instead  by  the  magnitude  of  the  Doppler-gradient  vec¬ 
tor,  or,  in  even  more  general  terms,  by  the  the  form  and  extent  of  the  wavenumber 
vector  or  k-space.  In  fact,  both  are  identical  only  for  the  very  special  case  of  side¬ 
looking  monostatic  SAR.  An  excellent  introduction  to  this  topic  can  be  found  in 
the  overview  article  by  Ender  [33].  Although  the  term  ‘Doppler  resolution’  might 
be  somewhat  misleading,  because  it  is  commonly  linked  to  Doppler  bandwidth,  we 
have  adopted  the  term  for  notational  convenience. 
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The  starting  point  to  determine  the  resolution  of  bistatic  SAR  is  the  three-dimensional 
geometry  depicted  in  Fig.  3.  The  effective  range  vector  from  the  transmitter  (from 
hereon  indicated  by  subscript  1 )  to  the  point  r  on  the  ground  and  back  to  the  receiver 
(subscript  2)  at  time  t  is  given  as  the  sum 

R(t,r)  =  5x(*.E)  +  Bz (t,L)  =  CSo.i  ~l)  +V.\t  +  (Bo,2~l)  +  v2t,  (60) 

where  R0  i  =  f?j(0,0)  are  the  ranges  from  the  scene  center  to  either  platform  at 
time  zero.  Determined  by  the  distance  between  the  transmitter  and  the  receiver,  the 
resulting  phase  history  of  the  received  signal  yields 

<Kt,r)  =  k(\\R1(t,r)\\  +  ||&(t,r)||),  (61) 

where  k  —  2n/X  denotes  the  wavenumber.  The  Doppler  frequency  is  given  by  the 
first  derivative  of  (61 )  as 

VD(t,r)  =  -  d4>^it  ^  = 

k _ v!i  {JM o,i  ~  L)  +  Hit) _ 

yj -  r)  +  UitY  (CSo.i  -t)  +  vj) 

+k _ ^((^o,2-Zl)+^0 _ 

\J {{Boa  -L)  +  Vat)'  ({Ro,2  ~  t)  +  JfcO 
=  k  (t,  r)  +  v2uE2  (t,  r))  ,  (62) 

where  =  a/\  |a|  |  is  the  unit  vector  in  direction  a.  For  simplicity,  the  unit  vector  in 
direction  R{  will  be  referred  to  as  Combining  the  scalar  k  and  the  unit  direction 
vectors  to  the  wavenumber  vector  fcf(f,r)  =  ku^t.r)  with  physical  unit  m-1,  the 
Doppler  frequency  yields 


Using  (63),  the  slope  of  the  Doppler  frequency  as  the  second  derivative  of  the  phase 

so(t,  r)  =  <p(t,  r)  is  given  as 

so(t,r)  =  +y2k2(t,r),  (64) 


which  can  be  further  written  as 

SD(t,r)  =  k^u^t^) +i/2u2(t,r)) 
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(65) 


,  /  (  £i  _ 

-1  KWMt^W  ||fi1(i,£)||3 

E2(t,r)R2(t,r)'v2 


+kv 2 


—2 


ILSa(t>£) 


ll£2(*,l)ll3 


%2PR2(^£)£2 

ll&(*.r)ll  +  ll&(«,r)ll  ’ 

where  describes  the  projection  operator  (matrix)  onto  the  complement  of  vector 

a. 

At  the  same  time,  the  gradient  of  <fi(t,r): 

V0(f,r) 


(to  dcj)  d<\> 
L^-’  ~dz~ 


in  (61),  at  point  r  can  be  calculated  as 


¥.<f>(t,r)  =  - 


kRi 


kR0 


ll&(MOII  ll&fct) 


=  -k^r)  -  k2(t,r). 


(66) 


The  gradient  of  the  phase  (distance)  function  determines  the  direction  and  the  max¬ 
imum  rate  of  change  of  distance.  Analogously,  the  gradient  of  0(f,r),  which  is 
identical  to  the  gradient  of  the  Doppler  frequency,  yields 


V^(f,r)  =  =  ~ki(t,r)  —  k2(t,r)  =  Vu>D(t,r),  (67) 


which  describes  the  direction  and  maximum  rate  of  change  of  Doppler  frequency 
for  any  point  r  on  the  ground.  Its  physical  unit  is  Hz/m.  From  (65),  the  first  deriva¬ 
tive  of  the  wavenumber  vector  vector  can  be  identified  as 


k(t,r)  = 


II  Ei(t,L) 


rp  i  (*>£)£<• 


(68) 


Note  that  the  Doppler  gradient  is  simply  given  as  the  negative  sum  of  the  two  in¬ 
dividual  wavenumber  vectors,  whereas  the  Doppler  frequency  represents  the  com¬ 
ponent  of  the  wavenumber  vectors  projected  onto  the  velocity  vectors.  The  same 
applies  to  the  Doppler  frequency  slopes  and  Doppler  gradient  slopes. 

4.2  Doppler-resolution 

From  this  k-space  analysis,  for  instance  [33],  it  is  known  that  the  resolution  (defined 
as  the  width  of  the  point  spread  function)  in  each  possible  direction  is  determined 
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by  the  reciprocal  of  the  maximum  extent/spread  of  the  k-set  in  this  particular  di¬ 
rection.  The  idea  of  using  the  gradient  to  calculate  this  maximum  extent,  i.e.,  to 
examine  the  bistatic  SAR  resolution,  was  originally  utilised  in  [34]  and  later  also 
in  [7,  35],  However,  these  authors  only  considered  the  wavenumber  vectors  at  the 
time  origin  and  neglected  their  change  over  time,  i.e.,  they  assumed  the  gradient  to 
be  constant  over  time.  This  is  a  rather  crude  approximation,  which  is  valid  only  for 
near-monostatic  (non-squinted)  configurations. 


Including  the  time  dependency  of  the  Doppler  gradient  (67),  the  maximum  k-set 
spread  in  Doppler  direction,  denoted  by  So,  can  be  calculated  as  the  norm  of  the 
vector  integral 


So(L)  =  - 


*ma: 

/ 


Vo;£)(f,  r)  dt 


(69) 


The  integral  must  be  computed  over  the  maximum  bistatic  observation  time  for  the 
considered  ground  point  r.  The  times  fmin  and  £max  denote  the  first  and  last  instant 
in  time  that  this  ground  point  is  ‘seen’  from  both  radar  antennas,  i.e.,  lies  simulta¬ 
neously  within  both  antenna  patterns.  For  radar  platforms  traveling  with  constant 
velocities,  the  bistatic  observation  time  was  calculated  analytically  in  Section  3. 


Using  (67),  this  becomes  in  closed  form 

^max 


Mr)  =  l 


J  +  kzitiL)  dt 


^  ||f^l  (^max>  Zl)  Ml(fmimll)  T  M2(^maxill)  M2(^min?  IL)  I]  •  (70) 


Hence,  the  exact  solution  for  the  resolution  in  Doppler  direction  AD  in  meter  is 
given  as  the  inverse  of  (70): 


A  M-A  _  C-l  _ 

A 

(1 1  ^ 

^DyL)  —  ~ 

(^maxt  L.)  IL]  (^min?  Il)  ll)  IL1?.  (^min?  Il) 

l/l) 

The  physical  interpretation  of  (7 1 )  is  that  the  Doppler  resolution  is  given  as  the  in¬ 
verse  of  the  sum  of  the  maximum  aspect  angle  variation  of  both  platforms. 


In  contrast,  the  maximum  Doppler  bandwidth  Bd  based  on  (63)  is  given  as 
k  k 

Bd(l)  71  Tyylkv  (ifl  (f  max>  Zl)  IL\  (^min  j  E))  “Fm  IL)  II))  > 

!|V]  ||  1  1 1  H.2 1 1 z  ~2 

(72) 


which  represents  the  sum  of  maximum  aspect  angle  variation  projected  onto  ci¬ 
ther  flight  direction.  In  the  monostatic  case,  with  Uj(tmax,r)  —  Ui(tmtn,r)  = 
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u2(t max,E)  —  U.2 (tmin.r),  which  we  will  call  a  for  convenience,  the  Doppler  band¬ 
width  is 

k  ollwll 

Bd(l)  =  -||w||uj,a  =  — F^-||fi||  3^2la-  (73) 

71  A  ~  ~ 

Equation  (73)  reveals  the  well-known  Doppler  bandwidth,  for  which  the  resolution 
is  calculated  through  A'D  =  ||v||/|£d|. 

From  (71)  it  is  evident  that  Bp  and  So  are  only  identical  if  =  1.  These 
two  vectors  are  only  collinear  for  a  special  symmetric  geometry,  namely  a  strictly 
side-looking  SAR. 


4.2.1  Special  cases 

To  gain  deeper  insight  into  the  ramifications  of  (71),  it  is  helpful  to  look  at  some 
special  cases  of  bistatic  geometry. 


4.2. 1.1  Monostatic  with  omnidirectional  antenna  pattern 


In  the  monostatic  case  when  the  two  direction  vectors  and  u2  are  identical,  i.e., 
u  =  u.i=  u2.  An  omnidirectional  antenna  diagram  results  in  a  theoretically  unlim¬ 
ited  observation  time  of  the  ground  point  r,  i.e.,  £min,max  =  Too  with 


,.  u  \  ,.  CBo.i  -r)  +  vt 

lim  u{t,r)  =  lim  — .  . -  ■  ...  = 

^±oo  t~±<x>  /i|^x  _r||2  +  2v,(R01  -  r)t  +  ||u||2 1 2 


=  ±T1\  ~ 


(74) 


As  expected,  the  unit  direction  vector  converges  to  the  unit  velocity  vector  for  infi¬ 
nite  time.  Hence,  inserting  (74)  in  (71)  leads  to 


Ad  = 


_ A _ _  A  _  A 

2||u(oo,  r)  —  u(— oo,r)||  2||2«„||  4’ 


(75) 


the  well  known  maximal  Doppler  resolution  of  a  quarter-wavelength,  independent 
of  the  scatterer  location. 


4.2. 1.2  Monostatic  stripmap  with  realistic  antenna  pattern 


The  maximum  observation  time  for  a  monostatic  geometry  is  given  as  the 
of  time  between  entering  and  leaving  the  mainbeam  during  the  overflight, 
times  are 


'"min, max 


=  =f- 


'3dB 


l&  -  dl 


period 

These 

(76) 
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which  in  turn  yields 


K:(^max?Zl)  min>2l)  — 


#3dB  V. 

2  Hull 


(77) 


V 


1  + 


(Bo 


Halt  IlSo-d 


^MB  +  ^ 


+ 


V1  + 


For  a  non-squinted  (or  slightly  squinted)  geometry  with  uRo_r  _L  v,  the  difference 
in  (77)  can  be  approximated  as  u(tmax,r )  —  M(£min,r)  ~  ^3do2^/ Hull which  results 
in  the  well  known  range-independent  (azimuth)  resolution 


Ad  = 


A 

l|2M3dB|| 


A 

2’ 


(78) 


where  /l  =  X/9-mb  denotes  the  real  aperture  of  the  monostatic  antenna. 


4.2. 1.3  Bistatic  with  omnidirectional  antenna  pattern 

Analogous  to  (74),  the  individual  direction  vectors  are  u,(o o,r)  —  «j(— oo,r)  = 
2 uv.,  which  inserted  into  (71)  leads  to 


2111^  +  uV2 1|  4\/l  4-  cosa 

where  a  denotes  the  angle  between  the  two  velocity  vectors.  It  is  apparent  that 
the  bistatic  Doppler  resolution  is  always  lower  than  that  of  the  monostatic  case, 
except  when  the  two  platform  fly  exactly  parallel.  In  contrast,  the  resolution  goes 
to  infinity  when  the  two  platforms  fly  anti-parallel,  i.e..  cosa  =  —1,  because  the 
aspect  angle  variation/spread  for  each  platform  cancels  out  exactly. 


4.2. 1.4  Bistatic  stripmap  with  realistic  antenna  pattern 


Let  us  assume  that  the  maximum  bistatic  observation  time  is  exclusively  determined 
by  one  of  the  platforms,  e.g.,  by  the  transmitter: 


to 


min, max  —  ^min,l,max,l  —  4~ 


#1  ||£o,i  L 


lliill 


(80) 


where  for  simplicity  of  notation  the  subscript  ‘3dB’  for  the  transmitter  beamwidth 
has  been  omitted.  This  assumption  is  not  restrictive.  In  contrast,  it  would  always 
be  fulfilled  in  cases  of  high-flying  (maybe  orbiting)  illuminators  and  relatively  low 
flying  receivers  with  ||i?0il  —  r II  ||^o.2  —  ill-  1°  order  to  keep  the  PRF  reasonably 
small,  i.e.,  to  avoid  range  ambiguities  while  maintaining  a  sufficient  swathwidth. 
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high  altitude  SARs  need  to  be  equipped  with  rather  large  antennas,  which  means 
Q\  -C  #2-  Analogous  to  (77),  the  direction  vector  difference  for  the  transmitter  is 
given  as 

'V 

Ml  (fmax,l  t  E)  Ml  (fmin.l  >  E)  ~  01  Tj  [7 1  (81) 

IK  II 


whereas  that  of  the  receiver  is 


M2(fmax,li  Zl)  M2(^min,liE)  — 

( 


(82) 


01 

H22o,i  —  e|| 

2 

IK  II 

—2 


,  IIM„,i-rll2fe||2e? 

1  ^  llfio,2-dlall2llli  4  / 


11^0,1  -Ell 
#1 - [1 — jj — M2, 

IK  II 


where  it  has  been  assumed  that  both  radars  work  in  quasi-stripmap  mode  with 
1.  Inserting  (81)  and  (82)  into  (70)  yields 


Bd  = 


ll^o,i 'Ell2, 


AIKII  V  11^0,2 -E 


—  r\\2  "—2 


+  2 


Ko,i  El 

11^0,2  -  El 


mK  +  IK  II2- 


(83) 


For  two  special  flight  geometries  this  term  can  be  further  simplified.  If  the  two  plat¬ 
forms  fly  exactly  perpendicular  to  each  other,  i.e.,  v\v2  =  0,  the  Doppler  resolution 
as  the  inverse  of  (83)  becomes 


(84) 


and  in  case  they  both  fly  parallel  or  anti-parallel,  i.e.,  v[v2  =  ±|KII  IK II  if  be‘ 


comes 


=  sign 


11—0,1  —  nil  IKII\ _ 

ILEo,2-eII  IK  11/  i± 


^4i 

llfio.i-dl  II112II  ’ 

!IMo.2-rll  fell 


(85) 


where  A\  =  A/<9i  denotes  the  real  aperture  of  the  transmitter.  In  (85)  it  is  interesting 
to  note  that  the  Doppler  resolution  will  be  better  than  that  of  a  monostatic  SAR  on 
the  transmitter  platform  when  the  ratio  between  distances  and  the  speeds  of  the  two 
platforms  are  chosen  appropriately,  e.g.,  \\Bo,i  -e||/IK II  >  II£o,2  —  eII/IKII-  The 
physical  reason  behind  this  is  that  the  combined  aspect  angle  variation/spread  of 
high  flying  illuminator  and  low-flying  receiver  can  be  larger  compared  to  that  of 
a  high-flying  monostatic  system.  In  [36],  this  result  was  derived  on  a  multi-page 
analysis  via  the  width  of  the  actual  point  spread  function. 


Fig.  12  shows  the  expected  bistatic  observation  time  and  the  corresponding  Doppler 
resolution  for  the  bistatic  airborne  scenario  listed  in  Table  6.  One  can  recognize 
that  the  bistatic  observation  time  is  non-symmetrically  distributed  along  the  ‘bistat- 
ically  illuminated’  area.  It  is  clear  that  the  Doppler  resolution  will  be  worse  for 
areas  where  the  observation  time  is  small,  such  at  the  edges. 
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Table  6:  Bistatic  airborne-airborne  configuration 


Parameter 

Transmitter 

Receiver 

Altitude 

Hx  =  7000  m 

H2  =  15, 000  m 

Velocity 

vx  =  130  k,.  m/s 

v2  =  220  m/s 

Course  angle 

I 

II 

o 

O 

$pi,2  =  45° 

Elevation  angle  to  scene 
centre  at  t  =  0 

$0,1  =  45° 

$0,2  —  35° 

Squint  angles  to  scene 
centre  at  t  =  0 

O 

II 

o 

O 

$0,2  =  90° 

Antenna  pattern 

sine 

sine 

Azimuth  beamwidth 

«3dB,l  =  3° 

^3dB,2  =  3° 

Elevation  beamwidth 

^3dB,l  =  20° 

^3dB,2  =  20° 

Wavelength 

A  =  0.0565  m 

Pulse  bandwidth 

B  =  50  MHz 

Sampling  frequency 

fs  =  50  MHz 

Pulse  Length 

T  =  5//s 

PRF 

680  Hz 

4.2. 1.5  Bistatic  stripmap  with  realistic  antenna  pattern  and  stationary 
receiver 

From  equations  (83)  to  (85),  it  is  evident  that  the  Doppler  resolution  for  a  stationary 
receiver  with  v2  —  0  is  identical  to  the  real  aperture  of  the  transmitter,  independent 
of  the  location  on  the  ground: 

&Z  =  A1,  (86) 

where  A  symbolises  the  moving  transmitter  and  the  stationary  receiver.  This  value 
is  exactly  one-half  of  the  maximum  value  for  a  conventional  monostatic  stripmap 
SAR.  Of  course,  the  bistatic  Doppler  resolution  can  always  be  enhanced  by  pointing 
the  antenna  beam  on  the  same  spot  on  the  ground  (spotmode)  during  the  flight,  and 
thereby  increasing  the  bistatic  observation  time,  and  in  turn,  the  extent  of  the  k- 
space.  In  contrast,  increasing  the  beamwidth  of  the  stationary  receiver  antenna 
will  not  enhance  the  Doppler  resolution  because  the  k-space  spread  is  exclusively 
determined  by  the  moving  transmitter.  However,  a  larger  receiver  beamwidth  will 
increase  the  imageable  area  on  the  ground. 

The  resolution  for  the  airborne-tower  scenario  introduced  in  Tables  4  and  5  is  deter¬ 
mined  by  the  transmitter  azimuthal  beamwidth  A£>  =  A/u3cjb.i  =  0.4516  m.  Since 
the  Doppler  resolution  for  a  stationary  receiver  is  constant  inside  the  bistatically 
illuminated  region,  it  is  not  explicitly  plotted. 

It  is  important  to  note  that  the  Doppler  gradient  (67)  is  only  varying  slightly  over 
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x-axis  [m) 

(a)  Bistatic  observation  time  in  seconds.  (b)  Doppler  resolution  in  meters. 

Figure  12:  Bistatic  observation  time  and  ground  Doppler  resolution  for  the  air¬ 
borne  transmitter  and  airborne  receiver  configuration  listed  in  Table  6. 


time  in  many  practical  cases,  such  as  side-looking  or  slightly  squinted  SAR  modes. 
In  such  cases,  k  can  be  treated  as  constant  over  time  and  can  be  taken  out  of  the 
integral  in  (70),  i.e., 

Bd  «  y  ll«i(0,l)  +lfcj(0,r)|| ,  (87) 

where  Tf,  =  tmax  —  tmin  is  the  bistatic  observation  time.  Replacing  0,  r)  with  the 
the  projections  (68)  at  the  time  origin  yields 


Bd*  ^ 


n 

A 


— - “Pd  (0,  r)v,  +  T7TT - -Pn 

Bo,i-d  ll^o,2 -dl 


PrMl)v2 


Tb 

A 


(I  -Mi(0»i)mi(0.e))^i 

llfiu  —  ill 


(I  -«2(0il)M2(0>l))l2 

11^0,2  —  ill 


(88) 


which  is  identical  to  the  result  originally  introduced  by  Cardillo  [34],  and  later  in 
[37],  For  only  slightly  squinted  geometries  where  u[(0,r)u„  *  0,  it  is  easy  to 
show  that  identical  results  as  those  for  the  previous  special  cases  (shown  above) 
can  be  achieved.  However,  (71)  represents  the  exact  solution  for  any  SAR  mode 
and  geometry,  and  should,  therefore,  be  preferred. 


4.2.2  Ground  plane 

In  most  cases  one  is  interested  in  the  resolution  on  the  ground  plane,  i.e.,  on  the 
separability  of  two  closely  located  points  on  Earth’s  surface.  In  order  to  accommo¬ 
date  this,  the  generally  three-dimensional  Doppler  gradient  vector  can,  for  instance, 
be  projected  onto  the  flat  ground  plane.  In  other  words,  the  maximum  change  of 
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Doppler  on  the  ground  is  along  the  projection  of  the  gradient  onto  the  ground.  Ap¬ 
plying  the  constant  projection  matrix 

/  1  °  °\ 

Pi  =  (I -s£)=  o  1  0  (89) 

\  0  0  0  / 

to  the  Doppler  gradient  in  (67) 

Vig)ujD(t,r )  =  Vc oD(t,r)  (90) 


and  following  the  analysis  form  (69)  to  (71),  yields  the  ground  resolution  of 


A  (®)  ( f  \  — 

A 

^Ml(^max?ZZ)  M.\  (^min  ?  ll)  1^2  max  ?  Il)  IL2  (^min  ■> 

.  1 ) 

4.3  Range-resolution 

The  idea  of  using  the  gradient  to  determine  the  bistatic  resolution  can  also  be  ap¬ 
plied  in  the  fast-time  direction.  Although  there  is  no  perspicuous  range  direction  in 
bistatic  geometry,  it  is  convenient  to  keep  using  this  term.  According  to  (1)  and  (3), 
the  total  travel  time  r  of  the  transmitted  pulse  to  a  scatterer  on  the  ground  (located 
at  r)  and  further  to  the  receiver  is  given  by 

r{t,r)  =  ^ (!!&(*, r) II  +  ||£2(t,r)||  ),  (92) 

where  c  denotes  the  speed  of  light.  The  gradient 

Vr(f ,  r)  =  -  («i {t,  r)  +  u2{t ,  r))  (93) 

c 

represents  the  direction  and  size  of  the  maximum  change  of  travel  time.  Although 
the  travel  time  is  strictly  speaking  dependent  on  time  (varying  position  of  the  radars), 
the  change  of  the  gradient  over  the  entire  aperture  time  is  typically  small,  so  that  it 
can  be  neglected  in  all  practical  cases.  Hence,  it  is  sufficient  to  analyse  the  gradient 
only  at  the  time  origin,  i.e.,  we  define  two  nearby  points  to  be  resolved  in  range 
as  long  as  they  are  separable  at  t  =  0.  For  a  monostatic  system,  for  instance,  this 
means  that  echoes  of  two  closely  spaced  scatterers  are  located  in  different  range 
cells  along  the  entire  synthetic  aperture  in  the  raw  data  matrix. 


Since  the  receiver  measures  different  arrival  times  proportional  to  the  system  band¬ 
width  B,  i.e.,  in  a  raster  of  spacing  St  =  1  /B,  the  range  resolution  can  be  written 
as 


Afl(r) 


1  IB 

II  Vt(0,  r)|| 


C 

B  ||t*i(0,r)  +M2(0,r)||’ 


(94) 
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where  Uj(0,r)  =  -  r)/IUSo.i  -  ill- 

Equating  the  two  directional  unit  vectors  in  (94)  confirms  the  well  known  constant 
slant-range  resolution  of  A n(r)  =  c/ (2 B)  for  a  monostatic  scenario.  In  the  bistatic 
case  (94)  can  be  re-written  as 


Aft(r)  =  —= - .  .  — 

\/2 B  \/\  +  cos  a(r) 


(95) 


where  a(r)  denotes  the  angle  between  the  range  vectors  from  the  point  r  on  the 
ground  to  the  receiver  and  transmitter,  respectively. 

4.3.1  Ground  plane 

Analogously  to  (90),  the  gradient  vector  (93)  must  be  projected  onto  the  two- 
dimensional  ground  plane  in  order  to  evaluate  the  ground  range  resolution  on  the 
flat  Earth’s  surface.  Therefore,  applying  the  projection  (89)  to  the  gradient  vector 
yields  the  general  form  of  the  bistatic  range  resolution 


Figure  1 3  shows  the  ground  range  resolution  over  varying  bistatic  angles  for  the  air¬ 
borne  transmitter  and  tower  receiver  geometry  used  previously,  see  Tables  4  and  5. 
The  loss  in  range  resolution  for  less  favourable  geometries  becomes  evident,  espe¬ 
cially  for  the  area  around  the  line  connecting  the  receiver  and  the  transmitter.  Here, 
the  resolution  decreases  drastically  from  roughly  30  cm  to  40  m.  However,  reason¬ 
able  imaging  (in  the  meter-range  resolution)  should  be  possible  for  pixels  further 
away  from  this  line,  even  for  a  forward-scattering  geometry  . 

4.3.2  Equivalent  bistatic  nadir  hole 

From  (96),  it  is  evident  that  the  ‘bistatic  nadir  hole’  is  given  at  the  location  on 
the  ground  r  at  which  the  norm  of  the  travel-time  gradient  vanishes.  Around  this 
location,  the  ground  range  resolution  tends  to  infinity,  i.e.,  meaningful  imaging  is 
impossible.  Re-writing  the  projection  matric  in  (89)  as  P.1  =  xsx'3  +  ygy's,  we  get 


In  order  for  the  norm  of  (97)  to  be  zero,  both  components  must  satisfy 


0 


0 


(98) 
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(a)  Quasi-monostatic  geometry,  bistatic  an-  (b)  Bistatic  angle  =  45°. 

gle  =  0  degrees. 
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(e)  Bistatic  angle  =  150°. 

Figure  13:  Ground  range  resolution  ii 
airborne  transmitter  with  vt 
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(0  Forward-looking  geometry,  bistatic  angle 
=  180°. 

meters  for  a  stationary  receiver  and  an 
ying  bistatic  angles  <J>0,2- 
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i.e.,  must  simultaneously  vanish  at  rg.  In  other  words,  the  vector  given  as  the  sum 
of  the  two  unit  vectors  pointing  from  the  location  ry,  to  either  radar  must  be  perpen¬ 
dicular  to  the  ground  plane,  i.e.,  Ug j(r0)  +  u0 2(zio)  J/.,- 

In  monostatic  SAR,  where  Mo,i(l)  =  Mo, 2(H)’  this  is  true  for  Lo  directly  beneath 
the  radar  platform,  which  is  typically  not  illuminated  by  the  antenna  pattern.  In 
bistatic  SAR,  however,  the  equivalent  nadir  may  shift  into  the  viewable/imaged 
area  depending  on  the  bistatic  geometry,  e.g.,  ||i?o a  »  i£,2||»  as  shown  in  Fig.  14. 
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5  Bistatic  SAR  raw  data  simulation 


To  develop  and  test  techniques  and  algorithms  to  process  bistatic  SAR  data,  it  is  es¬ 
sential  to  investigate  and  compare  their  fidelity  and  quality  on  simulated  data,  i.e., 
raw  data  sets  obtained  under  ideal  conditions.  The  developed  and  implemented 
point  target  simulator  comprises  all  necessary  components  to  be  realistic:  3D- 
platform  geometry  and  motion  vectors,  point  target  location  and  reflectivity,  imag¬ 
ing  geometry  (such  as  distances  between  illuminator-target-receiver  and  coordinate 
transformations)  as  well  as  antenna  pattern  and  variable  pointing. 

In  a  first  step,  an  airborne  bistatic  SAR  simulator  with  constant  velocities  and  al¬ 
titudes  and  a  flat  Earth  assumption  has  been  developed  and  coded  in  MATLAB. 
The  simulated  SAR  configurations  comprise  the  simple  stripmap  mode  as  well  as 
more  sophisticated  modes  where  the  antenna  is  continuously  steered  onto  the  scene 
center  during  the  flight  (spotmode,  sliding  stripmode,  etc.).  All  required  coordinate 
systems  and  angular  transformations  are  defined  in  appendices  D  and  E.  The  user 
inputs  are  the  basic  radar  system  parameters:  wavelength  A,  pulse  bandwidth  B, 
sampling  frequency  fs,  pulse  Length  T,  and  the  pulse  repetition  frequency  PRF,  as 
well  as  the  transmitter  and  receiver  parameters  at  the  time  origin  listed  in  Tables  4 
and  5. 


In  order  to  determine  the  maximum  raw  data  field  size  in  the  Doppler  direction, 
the  maximum  observation  time  T^ma*  is  calculated  via  (59)  for  each  pixel  within 
a  pre-specified  area  on  the  ground.  The  maximum  field  size  in  the  range  direction 
is  determined  by  the  longest  and  shortest  distance  to  any  pixel.  Half  a  chirp-length 
BT / 2  is  added  on  either  side  of  the  raw  data  field  to  prevent  focussing  errors  at 
the  edges  of  the  specified  regions.  The  pre-specified  region  on  the  ground  can  now 
be  sparsely  filled  with  P  scatterer  points,  where  each  point  is  characterised  by  its 
x,,ys-location  and  its  backscatter  amplitude  (radar  cross  section  (RCS)).  As  refer¬ 
ence,  the  bistatic  distance  to  the  center  of  the  scene  at  the  time  origin  is  calculated 

as  \m  —  iifluii  +  ii/m- 


The  actual  raw  data  generation  is  now  executed  in  a  straightforward  fashion,  namely, 
via  the  superposition  of  each  point  target  echo  at  each  pulse  transmission  time  (po¬ 
sition).  Instead  of  summing  the  reflected  echoes  in  the  time  domain,  it  is  com¬ 
putationally  more  efficiently  performed  in  the  Fourier-domain.  Here  the  time  de¬ 
lay,  which  requires  interpolation,  corresponds  to  a  multiplication  with  a  linear  fre¬ 
quency  ramp.  To  achieve  this,  the  ideal  transmit  pulse  in  baseband,  (2), 


s(ns)  =  exp 


for  ns  =  -JV./2, . . . ,  N./2,  (99) 


and  Ns  =  fi x(BT)  is  transformed  into  the  frequency  domain  S(ms)  =  P(s(ns)). 
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For  each  pulse  index  nv  the  range  vectors  f?j(np/PRF,r)  from  either  platform  to 
the  particular  scatterer  are  computed  via  ( 1 ).  Transforming  these  vectors  into  the  an¬ 
tenna  coordinate  system  with  (E.8),  i.e.,  by  application  of  Lap(<J>0,i,  ^o1i)Lp.s(^'p>t) 
and  using  equations  F.3,  yields  the  required  aspect  and  elevation  angles  to  this  scat¬ 
terer.  These  angles  are  then  used  to  calculate  the  antenna  pattern  weighting.  The 
delay  time  r(np,  r)  from  transmitter  to  the  scatterer  located  at  r  and  to  the  receiver 
is  given  as 


tM  =  ||&(VPRF,r)||  +  ||fi2(n,/PRF,r)||  -  ||fi,||  (|00) 

where  c  denotes  the  speed  of  light.  The  phase  multiplication 

S(ms)  •  exp  (  -  j27r/(ms)r(np)),  (101) 

where  f{ms )  denotes  a  frequency  bin  in  the  RF  domain,  includes  on  one  hand  the 
range  cell  migration  for  longer  synthetic  aperture  times,  and  on  the  other  hand,  also 
applies  the  Doppler  induced  phase  history.  Weighting  with  the  antenna  pattern, 
adding  up  the  returns  from  all  scatterers  and  applying  the  inverse  Fourier-transform 
to  the  result,  leads  to  the  desired  raw  data  signal  for  the  particular  pulse  index.  This 
procedure  must  be  repeated  for  all  pulses. 

Figure  1 5  (a)  shows  the  real  part  of  simulated  bistatic  SAR  raw  data  for  a  stationary 
receiver  scenario.  In  order  to  keep  the  raw  data  field  size  reasonably  small,  a  C-band 
transmitter  was  chosen;  the  system  parameters  in  Table  7  are  similar  to  those  of  the 
Environment  Canada  Convair-580  aircraft  [38].  Since  the  illumination  pattern  is 
not  symmetric  on  the  ground  for  this  perpendicular  configuration,  see  Fig.  1 1,  the 
slow  time  origin  does  not  occur  at  the  center  of  the  x-axis.  In  the  real  part  of 
the  raw  data,  typical  Fresnel-rings  are  visible  representing  ‘height  contours’  of  the 
two-dimensional  chirp  function.  Further,  the  ‘sinc’-antenna  pattern  weighting  in 
the  Doppler  direction  can  be  recognized  in  the  the  magnitude  plot.  Fig.  15  (b),  as 
well  as  the  range  migration  of  the  echo-pulse  through  various  range  cells. 

Currently,  a  spacebome  version  of  a  bistatic  SAR  simulator  is  under  development. 
It  is  based  on  a  fixed.  Earth-centered,  coordinate  system  and,  in  its  final  extension 
will  use  an  orbit  propagator  software  module  to  determine  the  satellite’s  state  vec¬ 
tors  at  each  pulse  transmission  time.  The  scatterers  are  placed  onto  the  curved  Earth 
surface  and  the  Earth’s  rotation  will  be  taken  into  account. 
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Table  7:  Simulation  parameters. 


Parameter 

Transmitter 

Receiver 

Altitude 

Hx  =  6900  m 

H2  =  20  m 

Velocity 

vx  =  132  m/s 

vx  =  0 

Elevation  angle  to  scene 
centre  at  t  =  0 

O 

LO 

TP 

II 

O 

$0,2  —  5° 

Squint  angle  to  scene 
centre  at  t  =  0 

0 

H 

0 

0 

0 

0 

0 

II 

<N 

o' 

Antenna  pattern 

sine 

sine 

Azimuth  beamwidth 

^3dB,l  =  3° 

^3dB,2  =  30° 

Elevation  beamwidth 

v3dB,l  =  20° 

v3dB,2  =  20° 

Wavelength 

A  =  0.0565  m 

Pulse  bandwidth 

B  =  50  MHz 

Sampling  frequency 

fs  =  50  MHz 

Pulse  Length 

T  —  5/xs 

PRF 

680  Hz 

(a)  Real  part.  (b)  Magnitude. 

Figure  1 5:  Real  part  and  magnitude  of  simulated  bistatic  raw  data  for  an  airborne 
transmitter  and  stationary  receiver  configuration,  \J>0)2  =  90°. 
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6  Bistatic  clutter  experiments 


As  part  of  the  TTRDP  (Trilateral  Technical  Research  &  Development  Program)  in¬ 
vestigation  of  bistatic  clutter,  the  MTI  experiments  team  conducted  a  bistatic  clut¬ 
ter  characterization  data  collection  campaign  in  the  summer  and  fall  of  2002.  The 
first  sets  of  experiments  included  the  so-called  Rooftop  test,  with  RADARSAT-1  as 
transmitter  and  the  AFRL  (US  Air  Force  Rome  Laboratory)  multi-channel  C-Band 
radar  as  stationary  receiver,  see  the  left  image  in  Fig.  16.  The  receiver  antenna  con¬ 
sists  of  five  separate  channels  that  allow  simultaneous  reception  of  different  signals 
in  parallel.  One  channel  is  designated  for  recording  the  direct  path  signal  from  the 
satellite  during  the  overpass.  The  right  image  in  Fig.  16  shows  part  of  the  planar 
antenna  array  as  well  as  the  additional  direct  path  receiver  horn  pointing  upwards. 
Since  this  horn  antenna  possesses  a  wide  beam,  it  does  not  require  an  adaptive 
azimuth  steering  capability  and  it  is  sufficient  to  tilt  it  to  a  constant  elevation  an¬ 
gle  according  to  the  flight  path  of  the  transmitter,  specifically  the  predicted  orbit  of 
RADARSAT-1 .  While  channel  one  records  the  reflected  signal  from  the  ground  that 


(a)  Rooftop-tower.  (b)  C-band  array  and  direct  path  re¬ 

ceiver  horn. 

Figure  16:  Stationary  bistatic  receiver  antenna  at  AFRL  in  Rome,  NY. 


is  received  over  the  entire  array  aperture,  channels  three-to-five  are  measuring  the 
echoes  over  different  subapertures  or  subarrays.  Channel  two  samples  and  records 
the  direct  path  signal.  The  array  is  rectangular,  3.5  m  in  width  and  1.2  m  in  height: 
each  of  its  eight  sub-arrays  is  0.45  m  wide. 
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6.1  RADARSAT-1  as  illuminator 


During  a  RADARSAT-1  bistatic  clutter  characterization  measurement  campaign  on 
January,  28,  2002,  several  Gbyte  of  data  were  collected.  A  detailed  data  description 
and  preliminary  data  analysis  were  reported  in  [39].  The  basic  system  parameters 
for  this  trial  are  listed  in  Table  8,  and  the  principle  bistatic  SAR  geometry  is  de¬ 
picted  in  Fig.  17  (not  to  scale).  While  the  satellite  is  flying  on  an  ascending  polar 

Table  8:  System  parameters  for  the  bistatic  clutter  trial  of  28  January  2002. 


Wavelength 

A  =  0.565  m 

Pulse  bandwidth 

B  =  11.56  MHz 

Pulse  Length 

T  =  42//s 

PRF 

1289  Hz 

Sampling  frequency 

fs  =  80  MHz 

Intermediate  frequency 

/IF  =  20  MHz 

orbit,  i.e.,  nearly  directly  towards  North,  it  is  illuminating  to  the  right.  The  station¬ 
ary  receiving  antenna  with  its  relatively  narrow  beamwidth  #3ciB  =  1°  is  pointing 
south-west  towards  the  town  of  Stockbridge. 


Figure  17 :  Bi static  geometry  for  the  trial  at  Rome,  NY  on  28  Jan.  2002.  The 


underlying  map  60  km  x  70  km)  was  generated  by  MapPoint 
© 2003  Microsoft  Corp.,  ©2003  NavTech  and  GDT,  Inc. 
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Since  the  Earth  is  rotating  eastwards2,  the  point  of  closest  approach  between  the 
antenna  and  the  satellite  tc  will  not  occur  at  broadside  4-  The  time  of  closest  ap¬ 
proach  is  at  the  apex  of  the  range  parabola,  as  illustrated  in  Fig.  17,  cf.  [40],  The 


(a)  Power  spectrum  of  received  IF  signal.  (b)  Power  spectrum  of  filtered  signal. 


(c)  Magnitude  of  received  signal.  (d)  Magnitude  of  filtered  signal. 

Figure  18:  SNR  enhancement  via  spectral  filtering  in  IF. 

recorded  data  stream  consists  of  the  sampled  real-valued  echo  signals  plus  a  super¬ 
imposed  1  MHz  calibration  tone.  Figure  18  (a)  shows  the  spectrum  of  part  of  the 
signal  after  rejection  of  the  calibration  tone  with  an  ideal  notch  filter.  The  roughly 
12  MHz  wide  transmit  pulses  can  clearly  be  recognized  around  the  Intermediate 
Frequency  (IF)  of  20  MHz.  Also  visible  are  discrete  interference  spikes.  While  the 
low-frequency  components  are  believed  to  be  caused  by  the  power  supply,  the  high- 
frequency  spikes  are  the  unwanted  higher  order  harmonics  of  the  calibration  tone. 

2The  velocity  at  the  Rooftop  receiver  latitude  is  approximately  340  m/s. 
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In  Fig.  18  (c),  a  small  portion  of  the  corresponding  time  signal  is  plotted  showing 
two  pulses  which  are  one  PRI  apart.  The  low  SNR,  typical  for  space-based  systems, 
is  apparent,  i.e.,  the  background  noise  is  relatively  high.  The  spectral  blanking  of 
everything  but  the  pulse  bandwidth  on  the  negative  side  of  the  spectrum,  see  Fig.  1 8 
(b),  results  in  a  significantly  increased  SNR:  the  pulses  are  completely  above  the 
background  noise  floor  in  Fig.  18  (d). 

After  spectral  shift  and  inverse  Fourier  transform  of  the  one-sided  spectrum  in 
Fig.  18  (b),  the  IQ-demodulated  signal  in  baseband  is  apparent,  with  Fig.  19  (a) 
showing  a  close-up  of  the  magnitude  and  phase  of  one  received  pulse.  The  nomi¬ 
nal  transmit  pulse  length  of  T  =  A2ps  and  the  Linear  Frequency  (LF)  modulation 
(parabolic  phase)  are  evident. 

The  data  were  then  range  compressed.  This  was  performed  by  convolution  with  the 
conjugate  complex  of  the  expected  theoretical  signal 

s(t)  =  exp  rect  ,  (102) 

using  the  nominal  values  for  bandwidth  B  and  pulse  length  T.  The  result  for  one 
particular  pulse  from  the  end  of  the  data  stream  is  plotted  in  blue  in  Fig.  19  (b). 
Although  the  resolution,  i.e.,  the  3dB  width  of  the  impulse  response,  is  close  to 


(a)  Magnitude  and  phase  of  pulse.  (b)  Compressed  pulse,  (blue)  theoretical  fil¬ 

ter,  (red)  matched  filter. 

Figure  19:  Compressed  pulse  sidelobe  reduction  via  matched  filter. 


the  expected  value,  the  first  sidelobes  are  very  high.  It  is  well  known  that  much 
better  performance  can  be  achieved  by  using  an  exact  matched  filter,  i.e.,  a  mea¬ 
sured  replica  of  the  transmitted  pulse  itself.  However,  since  such  a  replica  was  not 
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available,  one  of  the  first  pulses  was  extracted  and  used  to  compress  all  the  subse¬ 
quent  pulses.  The  red  curve  in  Fig.  19  (b)  represents  one  compressed  pulse  several 
seconds  (or  several  thousands  pulses)  later.  The  improvement  in  compression  is 
apparent,  particularly  since  the  first  sidelobes  approach  almost  perfectly  the  theo¬ 
retical  value  of  -13  dB  below  the  main  peak.  In  addition,  it  is  confirmed  that  the 
pulse  properties  are  stable  over  this  time  frame,  which  is  well  known  from  other 
experimental  SAR  systems. 

Figure  20  (a)  shows  the  magnitudes  and  the  Doppler  frequencies  (phase  differences) 
of  all  range-compressed  peaks  of  the  direct  path  signal  as  they  develop  over  time. 
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(a)  Magnitude  and  Doppler  frequency  along  (b)  Raw  data  array, 

time. 


(c)  Raw  data  array  after  range  migration  cor-  (d)  Range  and  Doppler  compressed  raw 
recti  on.  data. 


Figure  20:  Range  migration  correction  and  Doppler  processing  of  direct  path  signal. 
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The  blue  curve  is  the  amplitude  and  clearly  resembles  RADARSAT-l’s  antenna 
pattern  with  its  peak  at  around  7.5  s.  Strictly  speaking,  it  is  the  multiplication  of  the 
transmitter  and  receiver  diagrams;  However,  in  contrast  to  the  narrow  transmitter 
pattern,  the  receiver  pattern  is  very  wide  and  can  be  considered  as  almost  constant 
over  the  plotted  time  interval.  The  dent  in  the  centre  of  the  mainbeam  is  caused  by 
receiver  saturation.  However,  the  phase  information  is  not  effected  which  can  be 
seen  from  the  almost  perfectly  linear  Doppler  frequency  (red  curve)  in  this  region. 
The  Doppler  frequency  distribution  is  wrapped  at  the  sampling  limits  of  ±PRF/2. 
Large  phase  noise  (spikes)  occurs  at  areas  with  very  small  antenna  pattern  magni¬ 
tudes.  Re-ordering  the  continuous  data  stream  by  breaching  it  after  each  nominal 
PRI  leads  to  the  range-compressed  raw  data  array  in  Fig.  20  (b).  Due  to  the  chang¬ 
ing  distance  between  the  satellite  and  the  receiver,  a  significant  range  walk  can  be 
observed  for  the  direct  path  signal.  The  measured  peak  positions  for  each  pulse 
can  be  used  to  correct  for  this  walk,  which  is  commonly  called  range  migration 
correction.  Fig.  20  (c). 


In  order  to  compress  the  range-migration  corrected  raw  data  in  the  Doppler  direc¬ 
tion,  the  corresponding  Doppler  parameter,  i.e..  the  second  order  coefficient  in  (17), 
must  be  known,  or  can,  for  instance,  be  deduced  from  the  phase  history  in  Fig.  20 
(a).  It  is  either  possible  to  estimate  the  slope  of  the  Doppler  frequency  curve  or 
to  use  the  zero  crossings  to  compute  the  desired  Doppler  chirp  rate.  Using  the 
Doppler-time  relation  [41], 


ujD(t)  =  2 kvy  +  2k (t  -  4) , 

Kb 


(103) 


where  uy  denotes  the  cross-track  velocity  of  the  Earth,  and  inserting  the  times  trj  at 
which  the  zero  crossings  occur,  i.e.,  27r /d(4,;)  =  ±PRF/2,  we  can  solve  the  two 
equations  for 

PRF  RbX 

A  t  =  tr-U  = - (104) 

Using  the  time  differences  measured  from  Fig.  20  (a)  the  desired  acceleration  is 
estimated  as  v\l Rb  =  44.567  m/s2.  The  point  of  closest  approach  tc  is  charac¬ 
terized  by  the  time  at  which  the  Doppler  frequency  vanishes,  i.e.,  at  the  apex  of 
the  phase  parabola.  Setting  ujoitc)  —  0  in  (103),  the  across  track  velocity  can  be 
estimated  via  vy  =  v%/ Rb(tb  —  tc).  Since  the  direct  path  signal  is  under-sampled 
in  the  Doppler  direction,  multiple  ambiguities  arise;  Analysing  all  phase  maxima 
leads  to  two  likely  velocity  estimates  vy  =  299.376  m/s  and  vy  =  372.233  m/s, 
which  are  both  roughly  40  m/s  off  the  expected  Earth  rotation  velocity  of  340  m/s. 
This  discrepancy  is  caused  by  two  effects:  First,  the  satellite  is  not  in  an  exact  polar 
orbit,  but  has  a  slight  inclination  angle  so  that  the  radial  component  perpendicular 
to  the  pass  will  differ  somewhat;  and  secondly,  a  frequency  offset  between  the  two 
independent  local  oscillators  will  introduce  a  linear  phase  drift  as  shown  in  equa¬ 
tion  ( 1 5).  The  range  migration  corrected  direct  channel  data  can  now  be  focused  in 
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Doppler  direction  by  correlation  with  the  reference  function  (102)  using  the  deter¬ 
mined  chirp  rate.  The  result  is  plotted  in  Fig.  20  (d),  where  the  sharply  focussed 
2D-sinc  function  is  clearly  visible  together  with  one  ambiguity.  Since  no  weighting 
was  applied,  the  sidelobes  are  rather  high.  The  measured  3dB-mainbeam  width  is 
15  m,  which  corresponds  exactly  to  the  theoretically  predicted  value  in  (86)  (i.e., 
based  upon  RADARSAT-l’s  antenna  aperture). 

After  having  successfully  compressed  the  direct  channel  signal,  the  same  procedure 
was  applied  to  the  bistatic  channel.  Figure  21  shows  the  range  compressed  and 
range-migration  corrected  signals  from  channel  one,  i.e.,  the  bistatic  reflection  from 
Stockbridge  Hill  (as  measured  over  the  entire  antenna  aperture).  The  compression 
and  correction  was  done  with  the  information  gained  from  the  direct  path  signal. 
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Figure  21 :  Range  migration  corrected  bistatic  raw  data  field. 

Since  the  Rooftop  antenna  is  located  in  a  dip  ( 1 75  m  elevation)  it  is  actually  steered 
upwards  (with  an  elevation  beamwidth  of  2°)  towards  Stockbridge  (see  also  Fig.  24 
(b))  which  is  located  at  a  higher  altitude  (400m).  Due  to  this  geometry,  the  bistatic 
channel  also  includes  the  direct  path  signal  visible  as  the  strongest  reflection  at 
range  pixel  1020.  No  significant  difference  in  range  migration  between  near  and 
far  range  occurs  due  to  the  relatively  small  synthetic  aperture  time  of  about  one 
second.  Hence,  the  applied  correction  values  work  for  all  ranges;  the  trajectories 
are  straightened  out  everywhere  within  the  39  km  swath  (1500  pixels  xc/B).  The 
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roughly  45°  pointing  angle  of  the  receiver  can  clearly  be  seen.  Convolution  with 
the  reference  chirp  in  Doppler  direction  (102)  with  the  previously  determined  chirp 
rate  yields  the  final  image  Fig.  22. 

Again,  large  sidelobes  are  visible  because  no  weighting  function  was  applied.  Since 
the  receiver  beamwidth  is  relatively  small  (1°)  the  bistatic  image  (60  km  x  60  km 
size)  does  not  look  like  familiar  monostatic  images,  but  shows  only  a  thin  line. 
The  lack  of  details  is  simply  due  to  the  limited  ground  resolution  of  about  1 5  m  in 
Doppler  and  30  m  in  range.  A  close-up  look  of  size  2.2  km  x  7.8  km  in  the  vicinity 
of  the  receiver  location  is  plotted  in  Fig.  23  (a).  A  few  strong  scatterers  in  a  distance 
of  several  kilometers  are  evident,  which  most  likely  correspond  to  buildings  near  to 
AFRL. 


Time  [s] 

Figure  22 :  Range  and  Doppler  compressed  bistatic  SAR  data. 


Using  only  one  sub-aperture  of  the  entire  antenna  array.  Fig.  23  (b)  reduces  the 
antenna  beamwidth  to  about  8°  and  reveals  more  details  while  at  the  same  time 
losing  some  gain.  However,  as  expected,  this  will  not  increase  the  resolution  of  the 
system. 
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(a)  Entire  receiver  array.  (b)  Receiver  subarray  with  wider 

beamwidth. 

Figure  23:  Close-range  around  the  receiver  location. 

6.2  Convair  580  as  illuminator 

During  a  campaign  on  1  November  2002  over  the  AFRL  test  site  an  airborne  illu¬ 
minator  was  used  instead  of  RADARSAT-1.  The  radar  system  used  was  the  En¬ 
vironment  Canada  CV  580  SAR.  The  radar,  operating  at  C-band  (5.30  GHz),  was 
configured  to  use  a  transmitting  antenna  situated  on  the  belly  of  the  aircraft.  Hor¬ 
izontally  polarized  radiation  was  transmitted,  with  most  targets  of  interest  situated 
at  incidence  angles  of  approximately  47°.  The  Radar  system  was  run  at  a  PRF  of 
about  325  Hz.  However,  the  PRF  is  automatically  adjusted  to  match  the  aircraft 
velocity  during  the  flight  in  order  to  maintain  an  equidistant  ground  sampling.  At 
a  47°  incidence  angle,  and  an  altitude  of  approximately  6.5  km,  targets  of  interest 
were  offset  from  the  ground  projection  of  the  radar  flight  track  by  approximately  7 
km.  For  the  elevation  beamwidth  of  25°  the  swathwidth  on  the  ground  was  approx¬ 
imately  6.56  km.  Again,  Stockbridge  was  chosen  as  the  region  of  interest,  roughly 
28  km  east  of  the  Rooftop  receiving  antenna. 

Also  deployed  during  these  trials  was  an  active  bistatic  radar  calibrator  (ARC), 
Figure  24  (a),  which  was  located  on  Stockbridge  Hill  in  a  distance  of  30  km  from 
the  Rooftop.  The  two  horn  antennas  are  individually  steerable  into  different  di¬ 
rections;  while  the  receiver  horn  ‘looked’  towards  the  aircraft,  the  transmitter  horn 
was  directed  towards  the  receiver  tower.  The  ARC-recirculation  rate  was  set  to  1.8 
/is.  Figure  24  (b)  shows  the  surface  plot  of  a  Digital  Elevation  Model  (DEM)  for 
this  area  and  the  locations  of  the  Rooftop  antenna  (centre  cross),  the  ARC-reflector 
(cross  at  the  lower-left  of  centre)  and  the  flight  line  of  the  Convair  (yellow  arrow). 
The  DEM  illustrates  that  the  receiver  is  located  in  a  valley  looking  upwards  to  a 
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hilly  region.  The  flight  line  was  chosen  to  resemble  the  rooftop-RADARSAT-1  ex¬ 
periment,  i.e.,  a  quasi-forward  looking  bistatic  geometry.  These  trials  were  jointly 
coordinated  by  DRDC  Ottawa,  Canada  Centre  for  Remote  Sensing  (CCRS),  and 
AFRL  [42],  The  rooftop  receiver  radar  was  working  under  the  same  conditions 


(a)  Photograph.  (b)  Surface  plot  of  the  DEM  around  AFRL, 

NY,  (courtesy  R.  Hawkins  et.  al.  [42]). 

Figure  24:  Photograph  and  location  of  active  radar  calibrator  (ARC) 
at  the  1  Nov.  2002  airborne  trial,  crosses)  receiver  and 
ARC  locations,  arrow)  flight  path. 

and  with  identical  system  parameters  as  for  the  RADARSAT  trials.  The  airborne 
transmitter,  however,  emits  a  non-linear  frequency  modulated  pulse  of  T  =  7  ps 
duration  and  a  bandwidth  of  B  =  50  MHz.  However,  since  the  low-pass  filter  in 
the  final  stage  of  the  receiver  hardware  (before  sampling)  limits  the  signal  band¬ 
width  to  30  MHz,  a  resolution  loss  will  occur  compared  to  the  monostatic  case.  For 
these  parameters.  Fig.  25  shows  the  predicted  maximum  synthetic  aperture  times 
and  the  ground  range  resolution,  which  indicate  the  poor  range  separability  due  to 
the  forward  scattering  geometry.  Again,  one  of  the  first  measured  pulses  was  ex¬ 
tracted  in  order  to  compress  all  others  in  the  range  direction.  Fig.  26  (a).  The  pulse 
length  is  7  ps  long  as  expected.  The  magnitude  and  unwrapped  phases  history  of 
the  resulting  compressed  peaks  of  the  direct  path  is  plotted  in  Fig.  26  (b).  The  mag¬ 
nitude  is  the  result  of  the  multiplication  of  the  two  different  antenna  diagrams  of  the 
transmitter  and  receiver.  Due  to  the  particular  flight  geometry  the  direct  path  signal 
is  measured  in  sidelobe  regions  (both  in  azimuth  and  elevation)  of  the  receiver  and 
the  transmitter.  Therefore,  the  two-way  antenna  pattern  is  not  symmetric  and  shows 
high  sidelobes. 

To  re-order  the  continuous  data  stream  into  the  range-time  raw  data  array  the  vary¬ 
ing  PRF  has  to  be  taken  into  account.  The  range  migration  corrected  raw  data  field 
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(a)  Bistatic  observation  time.  (b)  Ground  range  resolution. 

Figure  25:  Predicted  bistatic  observation  time  and  ground  range  res¬ 
olution  for  the  Convair-Rooftop  trial  on  1  Nov.  2002 


is  shown  in  Fig.  28.  A  decibel  scale  was  chosen  in  order  to  emphasize  the  much 
weaker  ARC-returns.  Due  to  the  unfavourable  flight  geometry  (sidelobe  reception 
of  direct  path  signal)  any  range  migration  correction  using  the  direct  path  may  not 
be  reliable.  Instead,  it  was  found  that  it  is  much  more  precise  to  correct  for  range 
migration  by  using  the  ARC-returns  in  the  reflection-free  far  range  area  about  3 
km  distant  from  the  receiver,  because  the  ARC-signal  is  measured  within  the  main- 
beams  of  both  radars.  At  the  top  of  Fig.  28,  the  direct  path  signal  contained  in 
the  bistatic  reflections  can  also  be  seen.  The  relatively  strong  forward/sideward  re¬ 
flection  from  the  Hill’s  slope  appears  up  as  a  stripe  around  2  km  distant  from  the 
receiver.  The  bistatic  reflections  to  the  left  of  the  center  appear  weaker  than  on  the 
right  side,  which  is  in  perfect  agreement  with  the  DEM  structures  in  Fig.  24  (b).  The 
noise-shaped  reflection  in-between  are  the  reflections  from  the  relatively  flat  area 
between  the  Rooftop  and  Stockbridge.  Another  stripe  of  much  weaker  reflections 
is  seen  around  4  km  distance,  probably  stemming  from  a  second,  higher-located, 
ridge  behind  Stockbridge  (as  seen  from  the  receiver).  Multiple  re-circulations  from 
the  ARC  are  also  visible.  With  a  re-circulation  time  of  1.8  ps  and  a  sampling  fre¬ 
quency  of  fa  =  80  MHz,  the  number  of  samples  between  the  individual  signals 
should  be  nominally  144;  the  measured  value  is  145  pixels.  Each  pixel  corresponds 
to  a  distance  of  A r  =  cf  fs  —  3.75  m,  i.e.,  the  distance  between  two  circulations 
is  543.75  m.  Due  to  the  nearly  forward  looking  geometry,  the  entire  illuminated 
ground  swath  of  about  7  km  is  concentrated  (compressed)  in  a  roughly  500  m  wide 
stripe  in  the  ‘slant’  domain  (raw  data).  In  other  words,  the  ground  range  resolution 
is  drastically  reduced,  as  expected  from  the  theoretical  results  in  Section  4.3,  cf. 
also  Fig.  13  (e-f)  and  Fig.  25  (b). 
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(a)  Chirp  signal;  real  part  (red)  and  phase  (b)  Direct  path  signal;  magnitude  (red)  and 
(blue).  phase  (blue). 

Figure  26:  Received  chirp  and  direct  path  signals 
(magnitude  and  phase). 

In  order  to  compress  the  range-migration  corrected  data  in  the  Doppler  direction, 
the  Doppler  chirp  rate  must  be  known.  Unfortunately,  the  motion  data  for  the  Con- 
vair  aircraft  were  not  available.  Therefore,  an  autofocus  method  was  used,  based 
on  the  phase  information  of  the  ARC  in  the  far  range.  The  magnitude  and  phase 
history  of  the  ARC  at  3  km  range  is  plotted  in  Fig.  27  (a).  While  the  magnitude 
is  quasi  constant  over  a  large  region,  the  Doppler  phase  follows  a  typical  parabolic 
shape.  However,  the  apex  of  the  Doppler  phase  is  shifted  off  the  center  of  the  mag¬ 
nitude  pattern,  which  is  a  clear  indicator  for  the  presence  of  a  linear  phase  term, 
e.g.,  [40], 

The  technique  to  estimate  the  linear  and  quadratic  phase  coefficient  of  the  matched 
filter  is  adopted  from  the  moving  target  parameter  estimation  scheme  described  in 
[40].  There,  the  two  coefficients  are  influenced  by  the  along-track  (quadratic  term) 
and  across-track  (linear  term)  velocity  of  a  moving  object.  One  proposed  technique 
to  estimate  these  two  components  are  matched  filter  banks;  the  desired  parameters 
a  =  [ai,  02] 7  are  varied  throughout  a  bank  of  filters  of  form 

Stef  (t,  a)  =  exp  (j{ait  +  a2t2)) ,  (105) 

and  the  value  at  which  the  maximum  output  magnitude  occurs  represents  the  opti¬ 
mum  estimate.  Fig.  29  (a)  shows  the  filter  bank  outputs  for  the  chirp  rate  a2  varied 
between  50  s  2  and  200  s-2,  where  the  linear  term  «i  was  set  to  zero.  The  maxi¬ 
mum  appears  at  a  chirp  rate  of  121.4  Hz/s.  The  time  at  which  this  maximum  occurs 
is  tc  =  2.24  s,  which  represents  the  shifted  (incorrect)  position  of  the  ARC  in  the 
SAR  image.  Again,  this  displacement  is  due  the  presence  of  the  non-zero  linear 
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(a)  Magnitude  (green)  and  phase  (blue).  (b)  Magnitude  of  ARC-peaks. 

Figure  27:  Magnitude  and  phase  of  ARC-signal  at  3  km  range  and 
magnitudes  of  ARC-peaks  along  range. 


term  a\,  which  in  terms  of  GMTI  would  correspond  to  non-zero  radial  or  across- 
track  velocity  of  the  object.  As  known  from  the  GMTI  case,  around  the  optimum 
chirp  rate  the  matched  filter  responses  shift  away  from  tc  and  broaden;  the  ridge  is 
tilted  and  fans  out  at  the  edges  [40],  These  effects  are  also  caused  by  the  mismatch 
of  a!  in  the  matched  filter  and  the  data. 
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Figure  28  :  Range  migration  corrected  raw  data  in  logarithmic  scale. 
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Analogously,  a  matched  filter  bank  for  the  linear  term  component  ai  can  be  used. 
In  order  to  avoid  ambiguous  impulse  responses,  the  length  of  the  reference  function 
should  not  exceed  the  length  of  the  target  trajectory.  For  a  matched  filter  bank  map 
where  a\  was  varied  from  0  to  1000  Hz,  the  maximum  occurs  at  aii0pt  =  895  Hz. 
Since  the  ARC  was  stationary,  this  frequency  represents  the  difference  of  the  local 
oscillators  RF-frequencies  |/i  —  /2|  =  895Hz. 

Inserting  this  estimate  into  the  reference  function  and  repeating  the  a2  matched 
filter  bank  results  in  Fig.  29  (b).  The  ambiguity  is  resolved,  i.e.,  only  one  peak 
occurs;  the  ridge  is  now  correctly  shifted  to  the  broadside  time  4  =  -1.05  s  and 
the  fan  is  straightened-out  vertically.  The  impulse  responses  above  the  optimum 
value  (about  a2  =  117  Hz/s)  are  not  completely  symmetric,  which  is  caused  by 
the  presence  of  non-zero  cubic  phase  terms  [43],  This  term  is  most  likely  caused 
by  uncompensated  aircraft  motion  errors.  A  matched  filter  bank  for  the  cubic  term 
(not  shown  here)  resulted  in  a3,0pt  =  1.16  Hz/s-2.  Inserting  the  estimated  phase 
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(a)  Chirp  rate.  (b)  Linear  coefficient. 

Figure  29:  Matched  filter  banks  for  the  chirp  rate  and  the  linear 
phase  coefficient  based  on  the  ARC-signal  in  Fig.  21  (b). 

coefficients  into  the  reference  function  (105)  and  computing  the  convolution  results 
in  the  bistatic  SAR  image  shown  in  Fig.  30.  All  ARC-recirculations  are  correctly 
compressed  and  appear  as  a  “line  of  dots”  in  the  centre  of  the  SAR  image.  The 
linearly  decreasing  magnitude  (in  dB  scale)  along  range  is  confirmed  in  Fig.  27  (b). 

Due  to  the  very  small  ground  range  resolution  details,  such  as  streets,  buildings 
fields  etc.  in  the  Stockbridge  area  are  not  recognizable  However,  by  comparing 
Fig.  3 1  and  the  DEM  in  Fig.  24  (b),  certain  patterns  in  the  slope  structures  (different 
ridges)  in  the  DEM  show  up  as  reflections.  Further,  some  compression  gain  is 
observable  in  the  relatively  flat  region  between  the  receiver  and  Stockbridge  Hill. 
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Figure  30  :  Bistatic  SAR  image  for  the  1  Nov.  2002  airborne  trial. 
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Figure  31 :  Bistatic  SAR  image  for  the  1  Nov.  2002  airborne  trial  around 
ARC  location. 


-3-2-1  0  1  2  3 

Slow  time  [s] 


DRDC  Ottawa  TR  2004-190 


61 


6.3  Measured  phase  noise  characteristics 


This  section  investigates  the  phase  noise  characteristics  measured  during  the  bistatic 
clutter  experiments  at  AFRL  described  in  subsections  6. 1  and  6.2.  Since  the  direct 
path  signal  between  the  moving  transmitter  and  the  stationary  receiver  was  ob¬ 
served,  it  is  possible  to  explore  the  phase  noise  properties  under  practical  conditions 
and  to  verify  the  theoretical  results  introduced  in  Section  2.3. 

6.3.1  Rooftop-RADARSAT-1  experiment 

Figure  32  (a)  shows  the  direct-path  signal  phase  along  Doppler-direction  (blue 
curve).  This  phase,  confined  to  the  area  underneath  the  mainbeam  (T  =  6.9  s  to  8 
s),  is  given  as  the  integral  over  the  Doppler  frequency  plotted  in  Fig.  20  (a)  (red). 
The  green  and  the  red  curves  show  the  residual  phase  after  removing  a  second  and 
third  order  polynomial  fit  from  the  Doppler  phase.  As  expected,  the  satellite  shows 
only  small  motion  error  induced  phase  deviations  (e.g.,  a  relatively  small  cubic 
term)  over  the  synthetic  aperture  time  of  Ts  =  1.1  seconds.  The  IISLR  for  the  two 
curves  are  -15.3  dB  and  -19.3  dB  according  to  Table  2,  and  hence,  are  getting  close 
to  the  requirement  of  -20  dB.  Figure  32  (b)  confirms  this  finding  by  plotting  the 


(a)  Doppler  phase  (blue)  and  phase  noise  (b)  PSF,  ideal  (blue),  noisy  (red). 

(red,  green). 

Figure  32:  Measured  phase  noise  of  Rooftop-RADARSAT-1  experiment  and  im¬ 
pact  on  point  spread  function. 

point  spread  functions  with  and  without  the  shown  phase  noise  residual.  The  blue 
curve  represents  the  undisturbed  PSF  for  the  estimated  polynomial  coefficients,  i.e.. 
the  fitted  Doppler  rate.  The  red  curve  indicates  the  PSF  when  the  residual  noise 
(green)  is  added  to  the  signal  and  then  compressed  with  the  undisturbed  matched 
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filter.  Only  a  slight  one-sided  increase  of  the  peak  sidelobe  is  observed  due  to  the 
remaining  cubic  phase  error  but  virtually  no  increase  in  ISLR. 

6.3.2  Rooftop-Convair  580  experiment 

For  the  airborne  illuminator,  the  scenario  is  entirely  different.  Since  no  motion 
compensation  on  the  data  was  performed  due  to  non-availability  of  the  INS/GPS 
aircraft  data,  motion  induced  phase  errors  (caused  by  deviations  from  a  perfectly 
straight  flight  line  and  constant  velocity)  will  be  superimposed  on  the  oscillator 
phase  noise.  Figure  33  (a)  shows  the  residual  phase  noise  in  the  direct  path  signal 
over  eight  seconds  after  removal  of  quadratic  and  cubic  polynomial  fits  from  the 
measured  Doppler  phase  (blue  curve)  in  Fig.  26  (b).  A  much  larger  standard  de¬ 
viation  of  1.16  rad  and  0.75  rad,  respectively  is  seen  (these  values  are  too  large  to 
allow  computation  of  the  IISLR  via  (37)). 

The  drastic  impact  on  the  PSF,  even  for  the  third  order  phase  noise  residual,  is 
visible  in  Fig.  33  (b).  Alongside  a  significant  increase  of  the  peak  sidelobes,  a 
left-sided  increase  of  the  side  lobes  and  an  increased  ISLR  are  recognized.  While 
the  one-sided  sidelobe  increase  points  to  unresolved  higher-order  low-frequency 
phase  noise  terms,  the  two  peak  sidelobes  occur  as  ghost  targets,  which  might  be 
caused  by  the  sinusoidal,  low-frequency  structure  in  Fig.  32  (a),  [24],  However,  this 


time  [s]  time  (•] 


(a)  Phase  noise.  (b)  PSF,  ideal  (blue),  noisy  (red). 

Figure  33:  Measured  phase  noise  of  Rooftop-Convair  580  experiment  and  im¬ 
pact  on  point  spread  function. 

situation  improves  substantially  for  shorter  coherent  processing  times.  For  instance, 
halving  this  time  to  Ts  =  4  s  leads  to  much  smaller  standard  deviations,  i.e.,  to 
IISLRs  of  -14.2  dB  and  -19.7  dB,  respectively.  Correspondingly,  the  PSF  for  four 
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seconds  (third  order  polynomial  fit)  is  almost  ^distinguishable  from  the  blue  curve 
in  Fig.  33  (b)  and  is  therefore  not  plotted. 

Even  though  the  motion  error  induced  phase  noise  and  the  oscillator  phase  noise 
contributions  are  inseparable  without  proper  INS  data.  Fig.  34  reveals  some  deeper 
insight.  Here,  the  estimated  power  spectral  density  after  removal  of  up  to  second 
order  or  third  order  polynomial  fits  is  plotted  within  the  unambiguous  frequency 
band  (PRF  =  325  Hz).  The  Fourier  transform  reveals  discrete  phase  noise  at  about 
72  Hz,  which  are  known  to  be  caused  by  propeller  shaft  vibrations  of  the  Convair 
580.  The  other  spikes  at  about  30  Hz,  1 20  Hz  and  1 80  Hz  are  probably  due  to  frame 
vibrations  of  the  aircraft  (due  to  under-sampling  of  the  azimuthal  bandwidth,  some 
are  most  likely  ambiguous  frequencies  folded  back  into  the  PRF-region).  These 
vibrations  cause  the  oscillator  to  shake,  which  in  turn  will  perturb  the  oscillator 
frequency  (phase)  due  to  its  ‘g-sensitivity’.  A  detailed,  experiment  based,  analysis 
of  motion  and  vibration  induced  phase  jitter  and  the  resultant  limitations  on  the 
coherent  processing  time  will  be  performed  during  phase  II  of  this  TIF  project. 
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Figure  34  :  Power  density  spectrum  of  phase  noise  for  Rooftop-Convair  580 
experiment. 


In  contrast,  the  estimated  power  spectral  density  for  the  RADARSAT-1 -Rooftop 
experiment  shows  no  such  spikes,  as  expected. 
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7  Bistatic  SAR  processing 


Finding  fast  and  reliable  processing  algorithms  for  the  geometrically  more  demand¬ 
ing  bistatic  configuration  is  a  relatively  open  field  of  research  at  present.  Some 
work  on  this  field  had  been  reported  in  the  early  papers  by  Soumekh  [44,  45,  46] 
and  others  [47,  48],  However,  to  date,  most  bistatic  SAR  processors  are  based  on 
the  matched  filter  principle  and  are,  in  fact,  solely  time  domain  correlators  [49,  50]. 
Time  domain  correlation  or  convolution  is  universally  applicable  to  any  bistatic 
configuration  (including  monostatic)  and  is  relatively  simple  to  implement.  For 
these  reasons  we  have  adopted  the  matched  filter  processor  that  is  described  in  this 
section.  Future  work  will  certainly  include  investigation  of  alternative  processors. 

The  major  drawback  of  the  time  domain  technique  is  the  enormous  computational 
load  involved  that  leads  to  intolerable  processing  times.  Unlike  monostatic  SAR, 
the  kernel  function  for  bistatic  cases  is  non-separable  and,  hence,  the  use  of  fast 
operators  such  as  the  Fast  Fourier  Transform  (FFT)  is  prevented,  [16].  Recently, 
Ender  [5 1  ]  explored  the  theoretical  background  of  bistatic  Doppler  compression 
and  has  proposed  numerical  methods  to  separate  the  kernel.  Some  useful  analytic 
expressions  for  bistatic  geometry  can  also  be  found  in  [1 7] 

Let  us  assume  that  the  bistatic  geometrical  history  is  accurately  known,  for  instance, 
by  measurements  from  Global  Positioning  System  (GPS),  differential  GPS,  and 
Integrated  Inertial  System  (INS)  equipment  on-board  the  platforms.  Then,  the  range 
and  Doppler  history  of  any  scatterer  (location)  on  the  bistatically  observed  ground 
is  a-priori  known.  In  other  words,  the  range  trajectory  of  a  particular  scatterer 
can  be  predicted  within  the  two-dimensional  raw  data  field,  including  the  range 
walk  and  range  migration,  see  Fig.  35.  Note  that  the  range-trajectory  for  a  bistatic 
SAR  is  generally  not  a  parabola  as  compared  to  the  monostatic  case.  The  shaded 
areas  around  the  tracks  represent  the  duration  of  the  transmitted  pulses.  After  range 
compression,  which  is  done  identically  to  the  monostatic  case,  this  shaded  region 
is  compressed  down  to  a  narrow  line,  where  the  width  of  the  line  is  determined  by 
the  range-resolution  of  the  system,  see  Section  4.3. 

To  fully  identity  the  individual  scatterer  tracks  in  the  raw-data  array,  a  precise  time 
keeping  is  required  for  the  range  gating.  In  other  words,  one  must  be  able  to  assign 
the  exact  bistatic  distance  (near-range)  to  the  time  when  the  range  gates  are  opened 
at  the  receiver.  Since  illuminator  and  receiver  are  spatially  separated,  this  is  equiva¬ 
lent  to  a  non-trial  synchronisation  issue.  Time  synchronisation  can  actually  be  done 
directly,  i.e.,  with  an  extra  link  between  the  transmitter  and  receiver,  or  indirectly 
using  identical  stable  clocks.  The  first  technique  was  applied  during  the  bistatic 
clutter  trials  described  in  Chapter  6. 
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Figure  35  :  Illustration  of  a  bistatic  raw  data  field  (right)  showing  three  individ¬ 
ual  range  trajectories  for  scatterer  locations  specified  in  the  illumi¬ 
nated  footprint  on  the  ground  (left). 


Indirect  time  synchronisation  for  range-gating  based  on  the  pulse-per-second  (PPS) 
signal  of  the  GPS  offers  greater  flexibility  to  use  in  future  systems  as  well  as  reduced 
hardware  changes  on  the  existing  aircraft  [21,  28],  As  such  it  is  intended  for  the 
experimental  set-up  in  phase  II  of  this  project.  However,  the  fewer  modifications  to 
the  hardware  have  to  be  traded  with  less  accuracy,  i.e.,  a  larger  uncertainty  in  the 
range  cell  position;  this  error  is  expected  to  be  less  than  ten  ns. 

By  examining  Fig.  35,  the  task  of  how  to  compress  the  energy  for  each  image  pixel 
becomes  clear;  one  must  re-collect  the  spread  energy  along  the  a  priori  known  tra¬ 
jectories.  This  is  done  by  convolving  the  measured  data  along  the  trajectory  with 
the  theoretical  slow-time  signal  (matched  filter)  created  by  a  potential  scatterer  at  a 
particular  location  in  the  image.  A  slightly  different  way  of  executing  this  Doppler 
processing  is  referred  to  as  Algorithm  1 ,  as  described  in  appendix  G.  Here,  for  each 
range-line  (pulse  return),  a  complete  image  is  computed.  The  ensemble  of  images 
are  then  subsequently  added  coherently.  The  energy  of  the  echoes  of  each  pulse  arc 
distributed  throughout  the  image  according  to  the  a  priori  known  distance  between 
image  pixels  and  the  platform  positions.  The  introduced  algorithm  has  been  imple¬ 
mented  and  tested  on  the  experimental  bistatic  clutter  data  from  28  January  2002, 
as  described  in  Section  6. 1 . 

Figure  36  shows  the  Doppler  focussed  image  for  a  small  area  around  Stockbridge 
Hill.  Three  sperate  scattering  centers  along  a  quasi-straight  line  can  be  seen  along 
with  their  ghost  replicas.  In  order  to  understand  this  phenomenon,  three  equi- 
Doppler  (red)  and  equi-range  (yellow)  lines  are  superimposed  onto  the  image.  Ev¬ 
idently  the  ghost  images  appear  exactly  at  the  intersections  of  these  lines,  confirm¬ 
ing  ambiguous  echoes  caused  by  the  quasi  forward  looking  geometry.  Scatterers  at 
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these  locations  (intersection  points)  would  generate  exactly  identical  slow  time  his¬ 
tories  in  range  and  Doppler,  i.e.,  would  exactly  lay  on  top  of  each  other  in  Fig.  35 
and  possess  identical  Doppler  histories.  Further  details  of  this  algorithm  and  its 
properties  can  be  found  in  [52], 

In  contrast,  the  same  kind  of  ghost  images  did  not  appear  in  Fig.  23  since  a  direct 
convolution  with  the  estimated  Doppler  history  (chirp  signal)  was  applied,  which 
will  produce  only  one  single  match  whenever  the  matched  filter  lines  up  exactly 
with  the  measure  data.  However,  the  technique  used  in  Section  6  is  only  a  crude 
approximation  because  the  range  migration  correction  was  assumed  constant  over 
the  entire  image,  but  therefore  is  about  one-hundred  times  faster. 
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SAR  image  processed  with  algorithm  one  and  overt ayed  equi- 
range  and  equi-Doppler  lines  for  three  strong  scatterers. 


DRDC  Ottawa  TR  2004-190 


67 


8  Summary  and  outlook 


This  report  described  the  results  of  phase  I  of  the  TIF  project  on  bistatic  SAR. 
This  phase  lasted  one  year  and  was  concerned  with  gaining  a  deeper  knowledge  of 
bistatic  SAR  and  a  thorough  understanding  of  its  properties  based  on  theoretical 
analysis  as  well  as  analysis  of  existing  experimental  space-borne  and  airborne  SAR 
data  sets. 

A  complete  theoretical  framework  for  the  bistatic  SAR  resolution  on  the  ground 
was  developed.  Using  this  mathematical  tool,  the  range  and  Doppler  ground  res¬ 
olution  can  be  predicted  for  any  arbitrary  bistatic  flight  geometry.  By  using  the 
time-dependent  range  and  Doppler  gradients  and  the  bistatic  observation  time,  this 
can  be  achieved  without  actually  calculating  the  point  spread  or  impulse  response 
function,  which  makes  it  computable  in  quasi  real-time.  In  addition,  due  to  their 
simplicity,  this  set  of  equations  can  be  used  for  mission  planing  where  the  optimum 
flight  lines  must  be  determined  in  order  to  achieve  a  minimum  required  resolution 
for  a  given  area  of  interest.  Further,  performance  degradation  in  terms  of  reduced 
geometric  and  radiometric  resolution  of  the  point  spread  function  caused  by  severe 
phase  noise  due  to  separate  oscillators  was  investigated. 

In  order  to  verify  these  theoretical  results  a  complete  bistatic  raw  data  simulator  has 
been  developed  and  implemented  in  MATLAB®.  The  simulator  generates  the  raw 
data  with  high  accuracy  including  the  full  three-dimensional  geometry,  the  antenna 
diagrams,  and  the  complex  RCS  for  any  number  of  point  scatterers  located  on  the 
ground.  Hence,  realistic  effects  such  as  range  migration  and  magnitude  weight¬ 
ing  are  automatically  included.  A  more  practical  scene  including  features  such 
as  streets,  fields  etc.,  could  be  assembled  from  an  ensemble  of  point  scatterers. 
However,  the  computational  complexity  increases  drastically  in  proportion  to  the 
number  of  scatterers.  The  simulator  also  allows  antenna  squinting,  as  is  commonly 
done  in  special  SAR  modes  such  as  SlideSAR  and  SpotSAR. 

The  simulated  raw  data  are  used  to  test  the  performance  of  bistatic  SAR  processing 
algorithms.  In  particular,  a  ‘brute  force’,  time  domain  SAR  processor  was  intro¬ 
duced  in  Section  7.  Although  this  type  of  processor  is  very  time  consuming,  it  is 
based  on  the  matched  filter  approach  and  is,  therefore,  flexible  and  applicable  to 
any  complex  geometry. 

The  last  part  of  the  report  dealt  with  the  analysis  of  experimental  data  collected 
at  bistatic  clutter  trials  at  Rome,  NY  throughout  2002.  These  trials  comprised  a 
moving  transmitter  (either  the  RADARSAT-1  satellite  or  the  Environment  Canada 
Convair-580  aircraft)  and  a  stationary  receiver  antenna  located  on  the  premises  of 
AFRL.  Since  the  main  motivation  for  those  trials  was  to  collect  and  build-up  data 
bases  for  bistatic  clutter  properties  and  bistatic  GMTI  a  relatively  narrow  receiver 
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antenna  beamwidth  was  sufficient.  However,  due  to  the  relative  motion  between 
the  two  separate  radars,  these  data  sets  are  in  principle  also  suitable  for  Doppler 
processing,  i.e.,  bistatic  compression/focussing  in  the  ‘cross-range’  direction  to  en¬ 
hance  the  resolution.  Although  no  special  measures  were  taken  to  ensure  phase 
coherence  between  transmitter  and  receiver,  i.e.,  the  two  stable  local  oscillators 
were  running  freely  during  the  measurements,  both  data  sets  have  successfully  been 
focussed  to  the  theoretically  expected  resolutions.  This  was  mainly  due  to  the  rela¬ 
tively  small  synthetic  aperture  times  of  up  to  only  a  few  seconds  in  both  cases. 

In  Phase  II  of  the  project,  deeper  insight  into  the  functionality  of  a  bistatic  SAR  sys¬ 
tem  and  the  physical  properties  of  its  SAR  imagery  will  be  gained  by  experimental 
measurements  conducted  by  DRDC  Ottawa’s  airborne  SpotSAR  radar  system.  De¬ 
signing  and  constructing  a  secondary  receiver  which  will  be  placed  on  top  of  a 
stationary  platform  (e.g.,  a  tower  or  a  high  building)  and  locating  the  transmitter  on 
the  maneuvering  aircraft,  represents  the  simplest  useful  bistatic  configuration  and 
is  of  particular  interest  for  evaluating  the  performance  limits  for  scenarios  with¬ 
out  any  receiver  platform  motions.  Indirect  time  synchronisation  for  range-gating 
based  on  the  pulse-per-second  (PPS)  signal  of  the  global  positioning  system  (GPS) 
is  intended  for  the  experimental  set-up  because  of  the  reduced  hardware  changes  on 
the  existing  aircraft.  The  received  echoes  from  the  illuminated  ground  spot  would 
be  recorded  and  stored  for  off-line  processing.  Flight  trials  with  different  bistatic 
geometries  and  varying  imaged  scenarios,  such  as  different  vegetation  or  imported 
man-made  targets  (cars  etc.),  are  anticipated  to  be  conducted  during  phase  II  of  the 
project. 
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Annex  A 

Mean  and  variance  of  the  impulse  response 
function  in  presence  of  white  oscillator  jitter 


According  to  (32),  the  expectation  and  second  moment  function  of  the  impulse 
response  in  the  range  direction  are 


E  P(t) 

E|P(()|2 


,-JUlT 


T/2 


-T/2 
T/2  T/2 


^  J  J  gJ’rf  ((t-t'-tn-T)2-*'2)^^  ((t-t-tn-r)2-?) 
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•E  |ei(2i(t-*,-r)-S2(t-t')-Si(t-r-T)+Sa(t-i))J,  dt’di  (A  !) 


-T/2  -T/2 


Let  the  vector  E  be  Gaussian  distributed  with  expectation  zero  and  covariance  ma¬ 
trix  K,  i.e., 


—  —  [-i(fi))  — 2(^2))  “2(^4)]/  ~  -A/^OjK),  (A.2) 


where  the  covariance  matrix  is  given  as 


K  = 


0 

CSiSi^l 
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-  t3) 


°2 
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Ce2H2(^2  -  U) 


CE]Ei  (^1  -  £3) 
0 
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0  \ 

Ce2E2(^2  —  U) 

0 

°2  ) 


(A. 3) 


In  (A. 3)  the  independence  of  Ej  and  E2  has  been  utilised  as  well  as  the  property  that 
the  auto-covariance  is  an  even  function,  i.e.,  cee(£)  =  cee(~  t).  One  way  of  calcu¬ 
lating  the  expectations  in  (A.  1 )  and  (A.2)  is  through  application  of  the  characteristic 
function  < p(s)  of  E.  In  fact,  the  terms  can  be  expressed  as 


T>eUo)  =  Ee^-|  = 


2 


(A.4) 


For  the  mean  in  (A.l),  the  upper  left  2x2  matrix  of  K  and  s0  =  [1,  —1]'  must  be 
used;  for  the  mean  in  (A.2),  Sq  =  [1,  -1,  —1, 1]'  and  the  entire  K  must  be  used. 
Inserting  these  values  into  (A.4)  and  replacing  the  time  instances  by 


h  t  t  T ,  ^2  t  1  j  —~  i  t  T ,  /'4  t  t , 
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yields 


E  P(t)  =  e 


2,2 
a\  +a2 

2 


and  hence 

2  2 

|E  P(f)|  =  x i’7'  sinc(7r B(t  —  tn  —  T)). 


(A. 5) 


(A. 6) 


For  the  second  moment  function  we  get 

T/  2  T/2 

-T/2  -T/2 

(A. 7) 

Under  the  assumption  of  identical  oscillators  c^Ei  (i)  =  c=2=2(t )  =  cs=(0>  (A. 7) 
can  be  rewritten  as 


E|P(t)|2  =  e-W+^)^  |  J 


T/2  T/2 

E  |P(f)|2  =  e~2<j2  y  I  «-‘n-r)(t'-t)  ^ 

-T/2  -T/2 
T/2  T/2 

+e~2a2fr2  J  J  dt!dl  (A  8) 

-T/2  -T/2 

Providing  statistically  white  phase  noise,  i.e.,  che(0  =  cr2S(t),  where  6(t)  denotes 
the  Dirac  delta  function  ,  the  first  integrand  on  the  right  hand  side  of  (A. 8)  becomes 
the  line  integral 

2  T/ 2 

~T^~  /  (e2<T  ~1)v/2°^  (A.9) 

-T/2 

so  that 

E  |P(f)|2  =  ^e^2rj2  (e2*2  ~  l)  +  e~2a2  sinc2(7rB(f  -  tn  -  r)),  (A.10) 

and  hence,  with  (A. 5), 

a2p(t)  =  E  \P(t)\2  -  |E  P(t) |2  =  ^  (l  -  e~2°2)  .  (A.l  1) 
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Annex  B 

Standard  deviation  of  phase  noise  for  a  finite 
observation  time 


Let  E(t)  denote  a  wide  sense  stationary  stochastic  and  ergodic  phase  noise  process 
with  variance 


4  =  E  (=(()  -  E  =(f))2  =  E  =  (tf  -  (E  E(f ))2 ,  (B.  1 ) 

which  describes  the  variation  of  the  phase  around  its  mean  of  zero  for  an  infinite 
observation  time.  We  are  now  interested  in  the  variation  of  this  phase  process  (vari¬ 
ance)  for  finite  observation  times  T.  For  a  given  realisation  £( t )  of  the  stochastic 
process  measured  over  time  V  »  T,  we  can  estimate  this  variance  as  the  average 
of  the  differences  ((£r(£)  —  £r(f))2)  within  the  observation  interval  T.  Sliding  the 
observation  window  along  the  time  axis  we  can  see  how  this  (time-limited)  variance 
estimate  varies  over  time,  see  Fig.  B.l. 


T 


Figure  B.  1 :  Realisation  of  phase  noise  with  sliding  window  to  estimate  the  rms 
for  finite  observation  time  T. 


Mathematically,  this  can  be  expressed  as  the  convolution  of  the  phase  £(t)  and  the 
ideal  low-pass  filter  h(t)  =  1/T  rect (t/T),  i.e., 

<rf(f,T)  =  ^( t f  *  rect  (z{t)  *  rect  >  (B.2) 

where  *  represents  the  convolution  operator.  The  mean  value  of  this  variance  func¬ 
tion  over  the  entire  measurement  time  V  is 

T' /2 

°1{T)  =  ^  I  CT l(t,T)dt,  (B.3) 

-T'/2 

which  may  be  used  as  a  measure  of  the  rms  phase  noise  for  limited  observation 
time.  Applying  the  Fourier-transform  to  (B.2)  and  inserting  the  result  into  (B.3), 
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including  the  inverse  Fourier-transform,  yields 


T't 2  oo 


?(T)  =  J  j  (E(a>)  *  E(o>))  sine  (cuT/2)  exp( jcjf)  doudt  (B.4) 


-T'/2  -oo 
772  00 


—  —  J  J  (E(tu)  sine  (ouT/2)  )*  (E(u;)  sine  (ouT/2) )  exp(jout)  dou  dt, 


-T'/ 2  -00 


where  E(u>)  denotes  the  Fourier-transform  of  £(t)  (not  to  be  confused  with  the 
stochastic  process)  and  T  sinc(a;T /2)  is  the  corresponding  Fourier-transform  of  the 
rectangular  window.  Switching  the  integration  order  and  using  the  Dirac-function 
approximation  gives 

77  2 

^  J  exp  (j  out)  dt  =  sine  (ouT'/2)  =  S(ou)  for  V  1  (B.5) 

-T'/2 

and  writing  out  the  convolution  integral  leads  to 

OO  OO 

cr|(T)  =  J  J  E(ou  —  u)E(u)  dt/sinc  (ouT/2)  S(ou)  dou  (B.6) 

—  OO  — OO 
OO  OO 

- /  /  E(u>  —  u)  sine  ((cu  —  /2)  E(j/)  sine  (uT /2)  du  6(ou)  dou. 


— OO  — OO 


Exploiting  the  ‘sifting  property’  of  the  Dirac-function  and  the  fact  that  E(u>)  = 
E(— c u)*,  we  can  write 


<7|(^)  =  J  |H(w)|2  (l  —  sine  (ouT/2)2)  dou.  (B.7) 

—  OO 

In  order  to  investigate  the  statistical  properties  of  the  estimator  (B.7),  e.g.,  the  ex¬ 
pectation  value,  we  cannot  simply  replace  the  realisation  £(£)  with  the  correspond¬ 
ing  random  process  because  the  Fourier-transform  of  a  stationary  stochastic  process 
(which  has  infinite  energy)  does  not  exist.  However,  substituting  the  finite  Fourier 
transform  in  (B.4)  by  the  discrete  Fourier  transform 

1  A,_1 

£«  =  V  S  exp  U2nnk/N)  (B.8) 

k= 0 
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and  following  the  time-continuous  derivation  above,  leads  to 


1  N-l  N-l 

|H*,.|2  (l  —  sine2  ( nkK ))  =  ^hs(&)  (l  —  sine2  (irkK))  , 

fc=0  k=0 

(B.9) 

where  I^{k)  denotes  the  periodogram  of  length  N.  It  is  well  known  that  the  expec¬ 
tation  of  the  periodogram  for  sufficiently  large  N  is  identical  to  the  power  density 
spectrum  of  E,  i.e.,  E  /==(/:)  =  C—ik)  [30].  Finally,  in  time  continuous  notation, 
the  variance  can  be  computed  as 


<x|  =  j  Chs(^)  (l  —  sine2  (ujT/2))  du>  =  2  J  Che(^)  dio ,  (B.10) 

-oo  0.443 

T 


where  the  sinc-function  has  been  approximated  by  a  rectangular  window  with  a 
3dB-width  of0.443/7\ 


Alternatively,  the  same  result  can  be  calculated  as  the  expectation  of 


IfJ 

oo 

S(t)  -fj  E(s)ds 

2  1 

dt 

V  -oc 

—  OC 

J 

(B.l  1) 


where  both  integrals  of  stochastic  processes  are  defined  as  'mean  square'  integrals. 
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Annex  C 

Least  squares  estimation  of  phase  noise 
polynomial  coefficients 


In  contrast  to  (33),  let  E  be  a  discrete  stochastic  Gaussian  process,  i.e.,  En  = 
E(nAf)  for  n  =  1  where  AT  represents  the  sampling  period,  e.g.,  the 

PR1,  and  NAt  =  Ts  is  the  coherent  processing  time  (i.e.,  synthetic  aperture  time). 
Writing  this  discrete  stochastic  time  series  as  a  random  vector  E  =  [Ej, . . . ,  Ejv]', 
(40)  can  be  expressed  in  vector  notation  as 


E  =  A6  +  T,  (C.l) 


where 


A  = 

'  1  At  At2  At 3 

1  2A£  4A£2  8A  t3 

and  T  = 

'  Tx  ' 

T(A£) 

1  NAt  N2At2  N3At3 

.  T"  . 

T(NAt) 

(C.2) 


An  estimate  of  the  parameter  vector  b  =  [60, . . . ,  63]'  can  be  found  as  solution  of  the 
Least  Square  (LS)  minimisation  problem  min;,  |E  —  Ab\2  [29]: 

b  =  (A'A  )-1  A'E.  (C.3) 

Since  E  is  normal-distributed,  vector  b,  is  as  a  linear  transformation  of  E,  also 
normal-distributed  with  expectation  zero  and  covariance  matrix 

=  Ebb  =  (A'A)-1  A'E  EE' A  (A'A)-1  =  P#K=P#',  (C.4) 

where  P^  =  (A'A)  - 1  A'  and  the  covariance  matrix  of  E  is 


al 

chh(1) 

chh(2) 

•  •  •  cssLN) 

osd  i) 

aE 

chh(1) 

•••  cEE(A-l) 

cee(2) 

c==( 1) 

al 

•••  cm{N  -  2) 

.  (C.5) 

c-(A) 

ces(N  -  1) 

Czz(N  -  2) 

al 

The  factor  of  two  in  front  of  the  auto-covariance  function  (due  to  the  sum  of  two 
independent  oscillators,  see  (37)  has  been  dropped  for  clarity  reasons;  however,  it 
is  included  in  (41). 
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It  can  be  recognized  that,  for  instance,  the  first  column  of  K=  in  (C.5)  represents 
the  causal  part  of  the  phase  noise’s  auto-covariance  function  at  discrete  times  nAt, 
compare  (33). 

The  ith  element  of  the  parameter  vector  estimate  6*  =  is  also  normal-distributed 
with  zero  mean  and  variance 

=  “fe  =  u,P*KzP*Au„  (C.6) 

where  the  three-dimensional  unit  vector  ut  has  all  elements  zero  except  for  a  one  at 
position  i  =  1, 2  or  3. 

Subtracting  the  estimated  polynomial  from  the  original  phase  noise  vector  leads  to 

T  =  E-  Ab  =  E-  A(A'A)-1  A'H  =  P^H,  (C.7) 

where  =  I  -  A  (A' A)-1  A'  represents  the  projector  onto  the  complementary 
subspace  spanned  by  the  columns  of  A.  Again,  T  is  as  a  linear  transformation  of  a 
Gauss-vector,  normal  distributed  with  zero  mean  and  covariance  matrix 

Kt  =  E  XT  =  P^E  EE'Pi  =  PiK=Pi.  (C.8) 

The  variance  of  the  remaining  high-frequency  phase  noise  components  is  con¬ 
tained  in  the  diagonal  of  (C.8),  i.e.,  can  either  be  calculated  as  a\  =  1/A  tr{Kr  } 
( tr{  }  denotes  the  trace  operator)  or  as 

<4  =  y!nKrun  =  y/nPiK=Piun  =  a'„Ksan,  (C.9) 

where  an  =  P^Hn  represents  the  nth  column  of  P^  for  n  =  1, . . . ,  N. 

It  is  interesting  to  have  a  closer  look  at  (C.9).  The  covariance  matrix  K=  in  (C.5) 
belongs  to  a  important  class  of  Toeplitz  matrices,  the  so-called  circulant  matrices, 
where  each  column  is  a  ‘down-shifted’  version  of  its  predecessor  [53].  As  such, 
any  circulant  matrix  vector  product  constitutes  a  convolution  of  the  first  column  of 
the  matrix  with  the  vector,  which  in  fact  can  be  performed  via  three  discrete  Fourier 
transforms,  e.g., 

K=-an  =  c33(k)*an(k)  =  T~1{T{cE3(k)}  F{an(k)}}.  (C.10) 

The  ‘tilded’  quantities 

Css  =  [c==(l  :  A),  c=s(N  —  1  :  — 1  :  2)];  and  a„  =  [an,  0]'  (C.ll) 

in  (C.10)  represent  the  entire  (evenly-extended)  auto-covariance  function  and  the 
zero-padded  vector  a,  respectively  [53]. 
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Going  back  to  the  continuous  time  notation,  the  variance  in  (C.9)  may  therefore  be 
written  as 

OO  OO 

4  =  “n(t)-(css(t)*an(t))  =  I  an(t )  J  CssWAWeP**  <tfdt 

—  OO  —oo 

OO  OO  oo 

=  j  C~U)AU)  J  a Mel^'dt  if=  J  C~(f)\A(f)\2  dttC.12) 

—oo  —OO  —OO 

S - N/ - y 

A(-f)=A(f)- 

where  A(f)  =  ,F{d(£)}. 

Analogously,  the  variance  of  each  polynomial  coefficient  in  (C.6)  can  be  expressed 
as  a  filtered  version  of  the  spectral  density  function.  Let  the  elements  of  the  random 
vector 

Hi  =  Pa  Mi  =  A  (A' A)-1  6®  ^'  (C.13) 

represent  the  discrete  random  series  corresponding  to  the  time-continuous  stochas¬ 
tic  process  i>,(£).  Then 


OO 

4  =  J  C3s(f)  \Vi(f)\2df ,  (C.14) 

—  OO 

where  Vi(f)  is  the  Fourier-transform  of  Vi(t). 
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Annex  D 

Coordinate  system  conventions 


For  the  required  analysis  and  simulation  of  an  airborne  bistatic  SAR  system  it  is 
sufficient  to  concentrate  on  three  coordinate  systems  and  the  necessary  transforma¬ 
tions  between  them.  The  conventions  are  identical  to  the  ones  defined  by  Ender 
[54]. 

D.1  Scene  coordinate  system 

The  stationary  Cartesian  scene  system,  subscript  ‘s’,  describes  the  orientation  of 
the  underlying  terrain  and  serves  as  a  ‘master’  coordinate  system. 

xH  from  west  to  east  (left  to  right) 

y  cross-product  between  xs  and  zs 

za  positive  upwards 

The  origin  is  chosen  at  the  centre  of  the  scene. 

D.2  Platform  coordinate  system 

The  Cartesian  platform  system,  subscript  ‘p’,  describes  the  orientation  of  the  air¬ 
plane. 

xp  longitudinal  axis  or  roll-axis 
y  lateral  axis,  positive  to  the  right 
zp  cross-product  between  xp  and  yp,  positive  downwards 
The  origin  is  chosen  at  the  centre  of  mass  of  the  platform. 

D.3  Antenna  coordinate  system 

The  Cartesian  antenna  system,  subscript  ‘a’,  describes  the  orientation  of  the  antenna 
mounted  and  steered  on  the  platform. 

x 0  azimuthal  axis 

y  antenna  normal,  positive  to  the  right 


DRDC  Ottawa  TR  2004-190 


83 


elevation  axis,  cross-product  between  x a  and  y^,  positive  downwards 
The  origin  is  chosen  at  the  phase  centre  of  the  antenna  plane. 
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Annex  E 

Angular  transformation  conventions 


E.1  Euler  rotation  angles 


The  so-called  Euler  angles  are  defined  as  follows:  Initially  the  new  coordinate  sys¬ 
tem  is  identical  to  the  original  one.  The  following  sequence  of  operations  occurs: 

1.  Rotation  around  the  2-axis  through  the  azimuth  angle  'I'.  Counter  clockwise 
when  viewed  from  the  arrow  tip,  i.e.,  from  above.  The  two  systems  are  still 
aligned  in  the  2-direction. 

2.  Rotation  around  the  y-axis  (of  the  new  system)  through  the  elevation  angle  0. 
Counter  clockwise  when  viewed  from  the  arrow  tip,  i.e.,  from  the  left. 

3.  Rotation  around  the  x-axis  (of  the  new  system)  through  the  roll  angle  <3>.  clock¬ 
wise  when  viewed  from  the  arrow  tip,  i.e.,  from  in  front. 


1 1 


Figure  E.  1 :  Euler  rotation  angles. 


The  sequence  of  rotations  must  be  from  one  to  three.  The  individual  rotations  can 
be  described  through  matrices: 


(E.l) 
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and  the  composite  transformation  matrix  is  yielded  by 


L($,0,^)  =  L1($)L2(©)L3(^)  = 

(cos  ©  cos  $  cos  ©  sin  *P 

sin  <P  sin  ©  cos  'P  —  cos  $  sin  \P  sin  <P  sin  ©  sin  $  +  cos  (p  cos  ^ 
cos  <P  sin  0  cos  ^  +  sin  <P  sin  0  cos  <f>  sin  ©  sin  ^  —  sin  $  cos  'P 


(E.2) 


—  sin  0 
sin  $  cos  © 
cos  $  cos  © 


E.2  Angular  transformation 

The  new  coordinate  system  ‘N’  originates  from  the  old  one  ‘O’  by  rotation  around 

©NO,  'PnO  v*a 

LNo  =  L(<PNq,  ©no,  'Pno)  =  Li  (3>no)L2(0no)L3(^no)-  (E.3) 

An  arbitrary  vector  aQ  in  the  O-system  can  be  transformed  into  the  N-system  by 

On  =  LmoOo-  (E.4) 

The  inverse  of  Lno  is  identical  to  the  transpose  matrix,  i.e.,  L^o  =  L'no.  Let  ‘F’ 
be  a  further  coordinate  system.  It  yields 

Lpo  =  LfnLno-  (E.5) 


E.3  Scene  into  platform  coordinate  system 

The  angles  <Pps,  @ps,  \Pps  describe  the  orientation  of  the  platform  within  the  scene 
coordinate  system,  see  Fig.  E.2.  For  a  normal  flight  attitude,  the  angle  ©ps  is  chosen 
as  zero  and  <Pps  =  180°.  Hence,  the  transformation  matrix  reduces  to 

(cos 'Ppg  sin 'Pps  0  \ 

sintyps  cos  'P ps  0  ]  ,  (E.6) 

0  0-1  ) 

where  ^ps  denotes  the  course  angle.  Positive  ^ps  means  course  to  the  right  of  xs. 
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Figure  E.2  :  Angular  relation  between  scene  and  platform  coordinate  systems. 


E.4  Platform  into  antenna  coordinate  system 

The  angles  $ap,  ©ap,  'I'ap  describe  the  orientation  of  the  antenna  within  the  platform 
coordinate  system.  The  angle  ©ap  is  always  set  to  zero.  Positive  azimuth  angle  ^ap 
rotates  the  antenna  to  the  right  backwards  and  positive  depression  angle  <f>ap  inclines 
the  antenna  downwards,  i.e.,  first  rotating  then  tilting.  The  transformation  matrix 
then  yields: 

Lap  =  L($ap,0,®ap)  = 

(COS 'I'ap  sin  'I'ap  0  \ 

-  cos  $ap  sin 'I'ap  cos  'I'ap  cos  'I'ap  sin  $ap  ,  (E.7) 

sin  $ap  Sin 'I'ap  -  sin  $ap  cos  ^ap  eos<J>ap  / 


Figure  E.3  :  Angular  relation  between  platform  and  antenna  coordinate  systems. 
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Consequently,  the  transformation  between  the  scene  and  the  antenna  coordinate 
systems  is  given  as  the  product 

Las  —  Lap ( ^ap,  \hap)  '  LpS(^ps).  (E.8) 
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Annex  F 

State  vector  of  platforms  at  time  origin 


To  determine  the  location  of  both  platforms  over  the  ground  at  time  t  =  0,  it  is 
often  convenient  to  define  a  squint  angle  ^0,t  and  a  corresponding  elevation  angle 
$04  for  either  radar  antenna  i  =  1, 2  to  the  centre  of  the  scene.  The  antenna  normal 
vector  for  these  pointing  angles  is  determined  by  the  quiescent  pointing  vector  — 
yp  =  — =  —[0,  l,0]r  via  the  transformation  matrix  from  antenna  into  scene 
coordinates 

M&<  =  14 Uv  =  (LaP($o,i,  ^o, i)Lps(0))'ya.  (F.l) 

This  vector  expresses  the  unit  vector  from  the  scene  center  to  the  platform  at  the 
time  origin.  The  length  of  the  state  vector  can  now  be  determined  as  the  intersection 
point  of  with  the  x,  y-plane  elevated  to  a  given  altitude  Hi : 


—  Hi 


U  D 


Hi 


i±Q,i 


'LL  p 


(F.2) 


In  order  to  illuminate  the  centre  of  the  scene  at  t  =  0  with  exactly  the  peak  of 
the  antenna  mainbeam,  the  antenna  needs  to  be  steered  (squinted  and  tilted)  by  the 
directional  cosine  and  sine: 


Vo,i 


/  So,i 

-s\\RoA 

j  —0,i 

~S  II  /?  II 

II— 0,z  II 


(F.3) 


compare  (53). 
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Annex  G 

Bistatic  SAR  processing  algorithm 


Algorithm  1  Time  Domain  Processing  (courtesy  N.  Qotuk) 


1 :  for  /  =  1  to  L  do 

2:  for  n  =  1  to  TV  do 

3:  for  m  =  1  to  M  do 


4: 

5: 

6: 

7: 

8: 

9: 

10: 

11: 


{Loop  over  each  satellite  position} 
{Loop  over  image  cell  in  ^-direction} 
{Loop  over  image  cell  in  y-direction} 


Rsc  —  \/ ( (xs Xc)2  4“  (ys  Uc “h  (Zs  Zc)2 

_ {Calculate  satellite  to  cell  range} 

Rea  \J ((«Tc  Xq)2  “f*  {yc  Va)^  4"  (Zc  ~  Za)2 

{Calculate  cell  to  antenna  range} 

Rb  —  Rsc  “h  Rea 

{Calculate  bistatic  range  for  this  cell  and  satellite  location} 
Rsa  =  V((xs  Xa)2  “h  {ys  Va)2,  4“  {Zs  ~  Za)2 

{Calculate  direct  range  from  antenna  to  satellite  location} 
-^diff  =  Rb  Rsa 

{Calculate  range  difference} 

0  =  exp  (^i^diff) 

{Complex  phase  of  the  expected  return  echo} 

indexn  =  round  [^L] 

{Complex  phase  of  the  expected  return  echo} 

indexn  =  round 

{index  to  locate  the  range  cell  in  range  profile  (Ar  =  range  between 

each  sample  in  range  profile) } 

£nmlp  =  data(/, indexn)  •  conj(0) 

{Generate  a  set  of  n  x  m  pixels  for  each  satellite  position  by 
multiplying  complex  conjugate  of  calculated  phase  with  the  measured 

data} 

end  for 

{end  m  loop} 

14:  end  for 

{end  n  loop} 

15:  Enm  =  Enm  +  E^  p 

{Coherent  summation  of  images} 

16:  end  for 

{end  1  loop} 
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